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ABSTRACT 


Means for predicting the fluctuating pressures acting on 
externally blown flap surfaces are developed on the basis of 
generalizations derived from non-dimensionalized empirical data. 
Approaches for estimation of the fatigue lives of skin-stringer 
and honeycomb-core sandwich flap structures are derived from 
vibration response analyses and panel fatigue data. Approximate 
expressions for fluctuating pressures, structural response, and 
fatigue life are combined to reveal the important parametric de- 
pendences . 

The two-dimensional equations of motion of multi-element 
flap systems are derived in general form, so that they can be 
specialized readily for any particular system. An introduction 
is presented of an approach to characterizing the excitation 
pressures and structural responses which makes use of space- 
time spectral density concepts and promises to provide useful 
insights, as well as experimental and analytical savings. 
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EXCITATION, RESPONSE, AND FATIGUE LIFE ESTIMATION METHODS 
FOR THE STRUCTURAL DESIGN OF EXTERNALLY BLOWN FLAPS 



INTRODUCTION 


Short take-off and landing (STOL) aircraft concepts have been 
attracting much attention in the past several years, because 
of their potential operational advantages over more conventional 
aircraft. Of the several STOL configurations that have been 
given very serious consideration, those incorporating "externaily 
blown" flaps (e.g., see Fig. 1) have recently found increasing 
favor, and EBF aircraft technology currently Is the subject of 
extensive study and evaluation. 

Because the flaps of EBF aircraft are exposed to the direct 
impingement of the engine exhaust streams, as well as to the in- 
tense noise field that exists near the engines, the effects of 
the associated fluctuating pressures must be taken into account 
in the design of the flap structures. These effects are primarily 
of two kinds: (1) "sonic" fatigue of the flap airfoil structures, 

and (2) severe overall vibrations of the flap elements, with the 
associated high oscillatory loads in the structures (and actua- 
tors) that interconnect and support the flap elements. The pre- 
sent report is intended to provide some preliminary analytical 
approaches to assessing these effects, to revealing the import- 
ant parameters, and to suggesting improved approaches. 


The Sonic Fatigue Problem 

The complex problems of sonic fatigue life prediction and 
corresponding structural design fortunately may be simplified 
by considering them in terms of a sequence of sub-problems. 

These sub-problems consist of: (1) characterization of the fluc- 

tuating pressure excitation, (2) determination of the structural- 
responses to this excitation, (3) evaluation of the most signifi- 
cant associated oscillatory stresses, and (4) estimation of the 
corresponding fatigue life. 


FLUCTUATING PRESSURES ON EBF SURFACES 


Jet Configuration 

At several core-nozzle diameters aft of the nozzle exit 
plane, the flow field produced by fan-jet engines appears to 
be dominated by that due to the core jet (see Appendix A). 

It is reasonable therefore to estimate the fluctuating pres- 
sures produced by the exhaust from a fan-jet engine on the 
basis of the pressures associated with the core jet, and to 
make use of the extensive information available concerning 
ideal circular jets. 

The configuration of an ideal circular jet is sketched 
in Fig. 2. This shows a converging conical "potential flow 
region", surrounded by a diverging concial "mixing region". 

The total angle 2a subtended by the jet boundary typically 
is between 25 and 30 degrees (Ref. 3)* The length of the 
potential core is given (Ref. 2) by 

X =3.45 D ( 1 + 0.38 M) 2 , (1) 

c 

where D denotes the diameter of the (engine core) nozzle and 
M represents the Mach number of the exhaust stream. 


Velocity Fluctuations in Jets 


For estimation purposes it is also convenient to assume 
that the velocity fluctuations that are present in the jet in 
absence of an inserted flap are not altered substantially in 
the presence of the flap, and one may then interpret available 
velocity fluctuation data in terms of the pressure fluctuation 
information one requires. Figure 3 indicates how the axial 
turbulence intensity 



o 


I 


(2) 


Characterization of the fluctuating pressures on EBP sur- 
faces for design and estimation purposes is summarized in the 
first of the following sections and is discussed in some detail 
in Appendix A. The next of the major sections of this report 
deals with estimation of structural responses, of the associated 
stresses, and of the fatigue life of EBF structures (based on 
previously published information and on extensions of previously 
employed approaches), both for skin-stringer and for honeycomb- 
core sandwich structures. The assumptions underlying these es- 
timations, and their limitations, are also pointed out in these 
sections; a general , potentially extremely useful, approach toward 
characterizing the excitations and estimating responses is des- 
cribed in Appendix D. 


The Buffeting Problem 

This problem, which tends to be most significant at fre- 
quencies that are substantially lower than those of primary 
importance for sonic fatigue, also may be considered in terms 
of a sequence of sub-problems, namely: (1) characterization of 

the excitation, (2) determination of the associated vibratory 
responses, (3) evaluation of the most significant deflections 
and substructural loads, and (4) evaluation of these oscilla- 
tions, deflections, and loads in terms of ride quality and sub- 
structural design requirements. 

In contrast to the treatment of the sonic fatigue problem, 
this report deals with the buffeting problem in only a very 
preliminary manner. The final major section of this report 
presents the two-dimensional equations of motion of a general 
EBF system, in which the airfoil components are considered as 
rigid bodies interconnected by linear springs. These equations 
display the important parameters and provide a basis for carry- 
ing out natural-frequency and response calculations for specific 
EBF designs . 
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FATIGUE LIFE OF S KI N- ST RI NGER STRUCTURES 


Overview of Estimation Approach 

Conventional aircraft structures consist of skins, rein- 
forced by stringers, frames, and bulkheads (Fig. 7). Fluctu- 
ating pressures acting on the skins tend to induce complex 
vibratory deflections in the entire assembly, resulting in 
associated stresses, which — in turn — lead to structural 
fatigue . 

Because of the complexities of the excitations and re- 
sponses, currently available "sonic fatigue" design methods* 
are based on analyses developed on the basis of simplifying 
assumptions, coupled with empirically derived relations. These 
analyses in essence focus on one bay (i.e., one skin panel) at 
a time, ignore the complex spatial and temporal distribution 
of the exciting pressure by assuming the pressure always to be 
completely in phase over the entire panel, and compute the 
mean— square displacement response of the panel (mode by mode) 
to this spatially uniform, but time-wise random, pressure .They 
then calculate the maximum stresses from the panel modal dis- 
placements, and finally relate these calculated stresses to ex- 
perimentally measured stresses and fatigue data. 

The panel boundary conditions clearly play an important 
role; they not only affect the natural frequencies of the pan- 
el (which determine the parts of the excitation spectrum that 
dominate the response), but also the mode-shapes and therefore 
the relation between modal deflection and stress. Thus, pre- 
vious investigators have expended considerable effort on meth- 
ods for predicting the natural frequencies. 

In dealing with the panel responses and stresses, the rein 
forcing structures (i.e., stringers and frames) are considered 
essentially only as boundary conditions. They are in effect as 
sumed to deflect very little - an assumption that is likely to 


*Although these methods were developed to cope with the 
problem of fatigue induced by acoustic excitation, they may be 
expected also to be applicable (at least approximately) in many 
other cases of fluctuating-pressure excitation, including gener 
ally that due to impinging jets and tangential flows. In all 
cases, of course, the quality of the estimate depends on how 
well the actual situation matches the various underlying assump 
tions . 
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In view of Eqs. (7) and (17), then. 


L s = L 0A + 10 log 


2T • Hz 
1 + T 2 w 2 


( 22 ) 


Making use of Eq. (8), one may find the peak spectrum level. 


corresponding to 0 , to be given by 

ma x 


L s,peak ’ L 0A + 10 loE 


(*)- 


u 


10 log 


u 


+ c, 


ref/ 


with 


C = 10 log 


0.2 X 


«£=) ■ - - 


7 dB , 


(23) 


(24) 


where X re ^ represents a reference axial distance. To arrive at 
the above numerical value for C 13 X re f was taken as 10 ft, and 
the previously cited value of U r pf = 750 ft/sec was used. By 
combining Eqs. (20), (23), and (24) one finds 


L 


s ,peak 


/ U \ I y \ T +460 

152 + 30 log It]- 2 -) + 10 logU )~ 20 l0g — 520 (25) 

vref' '■ref' 


One similarly finds that the high-frequency spectrum level, 
corresponding to the high-frequency approximation <3> hi f req (f ) of 
Eq. (16), is given by 


J s,hi freq 


where 


L oa + 10 101 


fe )- 10 


c 2 = 10 log 


U _ • Hz 
ref 


,0.2tt 2 X_ -f : 


ref ref 


20 log 


rH ♦«. • 

r re f/ 


( 26 ) 


- 24 dB 


(27) 
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With this value of I, Eqs . (5) and (6) reduce to 

P 2 = (0.2*1 p q U 2 ) 2 ; (13) 

this then may be used to estimate the pressures on surfaces on 
which a jet impinges normally. 

From Figs. 5 and 6 one may similarly deduce that fgJL sur- 
faces on which jets impinge more nearly tangent ially , Vp 2 /q £ 0 . 12 , 
so that for conservative estimation purposes one may take 

p 2 ” = (0.12 q) 2 = ( 0.06 p Q U^) 2 . (1*0 

Clearly, the assumption of normal incidence leads to mean- 
square pressures that are higher by a factor of 16 than the 
pressures one obtains for more tangential incidence, and struc- 
tures that can withstand the normal incidence pressures for a 
given period may be expected in general to survive the tangen- 
tial incidence pressures for a longer period. 

Pressure spectrum (spectral density) . — The frequency 
spectral density * p (w) of the fluctuating pressures, for both 
the normal and tangential incidence cases, as has been stated, 
is given by Eqs. (7) and (8). The maximum value of the spec- 
tral density, which value is obtained for T^-u << 1, obeys 


$ (f) = 2p 2 T„ = ~~ p 2 * 0.0115 Xp 2 ^ , (15) 

max K f U o o o 


where the last expression has been obtained by substitution of 
Eq. (13). 

For high frequencies, on the other hand, — that is, for 
T f w >> 1, -Eqs. (7), (8), and (13) yield 

“F UP 5 ” P 

$ (f) « E = 2 = 0.029 -2-2. ( 16 ) 

hi freq 2ir 2 T f f 2 0.27r 2 Xf 2 Xf 2 

where f = oj/2tt denotes the cyclic frequency. 


7 


( 6 ) 
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1 
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PU 2 


represents the dynamic pressure at the exit. 

The frequency-spectral density $p(to) of the fluctuating 
pressure is shown in Appendix A to be of the same form as that 
of the fluctuating velocity component, and to obey* 


$ (to) 
P 


~ $ (f) 

2tt p 



T f /ir 
1 + T 2 co 2 


(7) 


where T„ represents a typical time scale (or inverse frequency) 
of the pressure or velocity fluctuations and obeys 

T * 0.1 X/U Q , (8) 


and where w denotes the radian frequency. 

From Appendix A one also finds that the pressure cross- 
correlation function <J> (s,x) for two points on the flap 

P 1 P 2 

surface near the jet axis, separated by a distance s, obeys 


Vp 2 (3 ’ t> ■ 


p- 


-s/L - T 
e e 1 


/Ti 


(9) 


where L denotes a length scale, called the correlation length, 
and is given by 


L » 0.025 X . 


(10) 


Near-Tangentially Impinging Jets. — For flap surfaces along 
which the engine exhaust flows essentially tangentially, the 
assumption of momentum flux annihilation would tend to overesti- 
mate the mean-square fluctuating pressure. From Figs. 5 and 6, 


*The spectral density ^(w) represents the mean-square 
pressure per rad/sec, whereas the spectral density $p ( f ) , ex- 
pressed in cyclic rather than radian frequency, represents the 
mean-square pressure per Hertz. 
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MAGNITUDE OF VELOCITY VECTOR, ft/sec 







NORMALIZED MEAN VELOCITY, U/U 0 


1.2 


t 1 1 1 1 n 1 r 



REDUCED RADIAL DISTANCE r/D 


FIG. 15 JET VELOCITY PROFILES (FROM REF. 27), COMPARED 
WITH BYPASS NOZZLE EXHAUST PROFILE. 
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C.G. OF 
n™ FLAP 



ATTACHMENT POINT BETWEEN 
nTH AND (n-l)TH FLAP 


CENTER OF ACTION OF 
AERODYNAMIC FORCES 
[COORDINATES: (X an ,Y an )] 


nTH FLAP IN 
EQUILIBRIUM POSITION 


INERTIAL REFERENCE 


COORDINATED SYSTEM ATTACHED TO TYPICAL 
FLAP ELEMENT, INDICATING LOCATIONS OF 
INTERCONNECTION AND AERODYNAMIC FORCE 
APPLICATION POINTS. 


FIG. 11 
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MACHINED 

LEADING 

EDGE 


NEUTRAL 
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SKIN 


TYPICAL HONE YC0M3- CORE SANDWICH FLAP ELEMENT. 



STRINGERS 
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FIG. 7 TYPICAL BAY OF SKIN-STRINGER STRUCTURE, AND TYPICAL STRINGE 
CROSS-SECTION, 
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VARIATION OF SURFACE PRESSURE FLUCTUATION ON A FLAT PLATE DUE TO JET 
IMPINGEMENT WITH LOCATION, AT X/D = 7 AND SEVERAL ANGLES OF INCLINA- 




VARIATION OF AXIAL TURBULENCE INTENSITY ALONG JET AXIS 
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corresponds to an intersection of the k b (w) curve, which is a 
plot of Eq. (6), with the w = online. Obviously, one may expect 
the product (integrand) to exhibit a very strong peak if the 
peaks of the factor functions coincide. The two frequencies at 
which this condition can occur are readily identifiable in the 
diagram: (1) at oo Q , for which k b = k Q , the admittance peak occurs 

at the same location as the jet noise peak. (2) At u)^, for which 

k = k, T , the admittance peak occurs at the same wavenumber as the 
b ri 5 

turbulent flow peak. 

Stresses. — Since the root-mean-square strain in a uniform 
beam or plate is very nearly equal to the ratio of the root-mean- 
square velocity to the longitudinal wave velocity in the material, 
calculation of the rms strain is a simple matter once one knows 
the mean-square velocity. Of course, one then merely needs to 
apply Hooke's law to obtain the rms stress. 

Methods are also available that permit one to account for 
the stress increases at boundaries, supports, or reinforcements. 
Determination of the motions of supporting structures (e.g., of 
plate ribs or frames) and of the associated stresses involves 
additional calculations, which often can be formulated in rela- 
tively simple terms. 
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The (space and time) Fourier transform of Eq. (4), when 
written in terms of the velocity v instead of the displacement 
y, may be shown to be 

(k 4 - k 4 ) v = p(k,w) , (5) 

where k and co represent the wavenumber and radian frequency 
(which replace x and t as a result of the Fourier transform 
process) and p(k,oo) represents the Fourier transform of p(x,t) 
— i.e., the spectral pressure amplitude. The wavenumber k^ 

obeys 

k 4 = w 2 m/B (6) 

b 

and is that wavenumber at which free bending waves travel with 
out dimunition (in the absence of damping). 

The admittance Y is defined as the ratio of the spectral 
velocity to the spectral pressure; for the beam under consider 
ation here one finds 


(7) 


As evident from this equation or from the sketch of this func- 
tion appearing in Fig. 21, the admittance has a sharp peak at 
k = k^. In fact, this peak rises to infinity for undamped 
systems, but remains finite for realistic structures that al- 
ways have some damping. 

Admittances of Other Structures . — Equation (7) probably 
represents the simplest admittance function of practical in- 
terest. However, expressions are also available for the ad- 
mittances of beams with uniformly spaced masses attached, or 
of one-dimensional plates with uniformly spaced ribs. These 
expressions are more complicated than that for a uniform beam, 
but are no different in concept. 

For two-dimensional structures, the admittance functions 
involve two wavenumbers, corresponding to Fourier transforma- 
tions on the two spatial coordinates, but otherwise again do 
not differ in concept from the simple beam admittance. 


Y (k , id) = 


v 


com 

ii - !k)' 


' k 4/ 

n 


7 ^ 


$ ( k , co ) = | p ( k , co ) | 2 L“ 1 T~ 1 


( 3 ) 


where L and T, respectively, denote the spatial interval (here, 
in one dimension for this introductory one-dimensional discus- 
sion) and the time interval over which the pressures are sampled. 

Surface Pressures in Turbulent Flow . - Figure 19 indicates 
schematically the typical behavior of the pressure spectral den- 
sity $(k,io) associated with turbulent flow along a surface. If 
one considers the curves of constant spectral density as "contour 
lines" in the k,w plane, one notes that the $ hill has a ridge 
along a line whose slope is co/k = Uoo, where Uoo represents the 
speed of the flow along the surface. For blown flap surfaces, 

Uoo generally is considerably less than the speed of sound c n , as 
also indicated in the figure. 

All points on the w = kUoo line represent energy travelling 
at the speed of the flow. If all of the energy would travel at 
the flow speed, — i.e., if the turbulence were "frozen" into 
the flow and would convect at the flow speed — the contour lines 
would collapse into the Uoo line. The spread of the contour 
lines about this line reflects the "unfrozen" nature of boundary 
flow, — i.e. , the presence of a distribution of components tra- 
velling with different speeds. 

At any particular frequency Uj, the component that travels 
with the speed of the flow has a certain wavenumber kj.j . This is 
called the hydrodynamic wavenumber at that frequency, and it is 
the wavenumber at which there exists the greatest fluctuating- 
pressure energy at that frequency. Similarly, the component at 
Wj that travels with the speed of sound has a wavenumber k 0 , the 
acoustic wavenumber. Sound waves travelling (at grazing inci- 
dence) along the surface correspond to the c Q line of the figure. 

The lower part of Fig. 19 shows a "slice" taken through the 
upper plot at the constant frequency a^. The peak in the spec- 
tral density at the hydrodynamic wavenumber k H again displays 
the fact that most of the energy travels at the speed of the 
flow. It is important to observe, however, that significant 
amounts of energy also travel at other speeds, particularly at 
sDeeds near Uoo. Indeed, the fact that $(<i)j,k) is finite for a 
range of values of k below k Q indicates that some energy travels 
at sceeds greater than the sound speed. 
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first of these functions is known to acousticians (and to radar 
and sonar engineers) as the directivity pattern of the source 
of excitation; it represents the spectral density of the sur- 
face pressures. The function representing the structural stiff- 
ness is known to structural dynamicists as the surface (spectral) 
admittance of the structure, and to acousticians as its array re- 
ceiving pattern. The available physical understanding of these 
functions should facilitate direction and evaluation of research 
and design programs. Since only integrals of products of func- 
tions are required, only minimal computational difficulty is in- 
volved. 

Experimental Savings and Accuracy Gain . — With the STSD ap- 
proach, the (spatial) spectral density of fluctuating pressures 
acting on a structural surface may be measured by means of an 
array of pressure sensors flush-mounted on the surface; the out- 
puts of ail of the (spatially distributed) sensors are sampled 
simultaneously and recorded. The records are later digitized 
and entered into a digital computer for calculation of the spec- 
tral densities. Only one experimental run is needed for each set 
of independent variables. 

Alternate approaches to representing surface loads in es- 
sence require the repositioning of pairs of transducers and many 
reoetitive runs for each set of independent variables. The re- 
corded sensor outputs from all such runs typically are cross- 
correlated for a range of time delays, leading to space-time 
cross-correlation functions. Because data from different runs 
are obtained and analyzed separately in these approaches, loss of 
accuracy results due to random experimental errors. The STSD 
approach, on the other hand, is based on a single run — not 
only reducing this source of errors drastically, but also lead- 
ing to a considerable saving in run time. 

Extensions . - The STSD approach is ideally suited for 
dealing with the responses of surface (e.g., skin) structures 
to snatially homogeneous excitation. Where the responses of 
interior or supporting structures (e.g., ribs, stringers) are 
of interest, these may be expected to be determined relatively 
simDly from the surface structural responses. 

As the spatial inhomogeneities of the exciting pressure 
field and of the structure increase, the STSD descriptions re- 
quired to provide sufficient accuracy become more complex. In 
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TABLE C-I. - FATIGUE LIFE CORRECTION FACTORS FOR ALUMINUM ALLOYS 
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■^Completely reversed flexural stress that flat specimens can endure for 10 s cycles. 
From Table 3-3*l(c), Ref. 34. 
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take the endurance limit as the fatigue stress. Aluminum alloys 
do not have endurance limits, in general (Refs. 34 and 35); for 
such alloys one needs to define S as corresponding to any fixed 
number of cycles, say N = 10 8 . 

Correction for Fa tigue Stress . - If one summarizes Eqs . (45), 

(52), (71) and (75), "as 


N = B(a l /a ref ) 6 , 


then one may write 


ref 


h:-(S) 


-i/e 


and, introducing Eq . (1), 


S 2 a j 


a ref S i °ref i 


S 2 /n\ 

S \b/ 


-i/e 


( 2 ) 


(3) 


(4) 


Then , 
found 


the number N of 
to obey 


N = B 


cycles that material 2 can withstand 



is 


(5) 


where Ni denotes the number of cycles one calculates for the 
basic material — i.e., by use of Eq . (2) — and 


k = (S,/S ) B 
m 2 1 


( 6 ) 


is a correction factor that accounts for differences in the 
material’s fatigue properties. 

Since the fatigue life is proportional to the number of 
cycles to failure, one may obtain the fatigue life of a struc- 
tural component of any aluminum alloy by multiplying the life 
one calculates for that component on the basis of Eqs. (45), 
( 52 ), ( 71 ), (75), by the appropriate correction factor k m , as 
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where E represents the modulus of elasticity of the skin material. 
The shear stress In the core obeys 


* - - - 5 v - ! (c) 3 ». ) • (11) 

Thus, the maximum shear stress t , which may be seen to occur at 

3 ITla.X 

the ends of the beam, is related to the maximum skin stress a max 
due to flexure (which occurs at the beam center) as 


T max 

a 

max 


ttB 

LAcE ' 


( 12 ) 


Rectangular Section Beam . - For a full-depth honeycomb beam 
with a rectangular cross-section of width e and thickness H, and 
with facing sheets (skin) of thickness t s , 

B = EH 2 et /2 , A ~ eH . (13) 

s 

Then the correction factor term appearing in Eq . (8) becomes 


n 2 B _ 7 t 2 E ^ s 

L 2 AG G L 2 


Since Ht s /L 2 generally is very small, the expression of Eq. (14) 
may be exoected to be small compared to unity, except for very 
soft core's, for which tt 2 E/G is very large. Thus, except for such 
soft cores, the bending natural frequency expression of Eq. (9) 
may be expected to suffice, since the factor by which is multi- 
plied in Eq. (8) then is very nearly equal to unity. 

By substitution of Eqs . (13) into (12) and noting that 
c = H/2 one finds that for a rectangular section beam 


TmaX * Is (15) 

a L 

max 


6b 


APPENDIX B 

SHEAR EFFECTS IN HONE YCOMB- CORE SANDWICH BEAMS 


Since the honeycomb core, when considered as a continuum, has 
a relatively low shear modulus, shear deflection may play an im- 
portant role in the dynamics and stresses of honeycomb core sand- 
wich beams. This appendix presents an approximate analysis of 
the effects of the honeycomb beam shear stiffness. 

In Ref. 26 there is summarized an analysis of the stress re- 
sponse of simply supported honeycomb-core sandwich panels. The 
analysis presented below proceeds in the same manner as that of 
Ref. 26, but applies for beams instead of panels. 

Relation Between Flexural and Shear Deflections . - The total 
deflection w of a beam may be considered as composed of a compon- 
ent w^ due to bending and of a component w s due to shear, that is, 

w = w, + w ( 1 ) 

b s 


For a beam vibrating at its fundamental resonance frequency w, 
one finds from simple beam theory that the deflections must sat- 
isfy the differential equation 


where the primes indicate differentiation with respect to the 
longitudinal coordinate x, B denotes the flexural stiffness of 
the beam, and y its mass per unit length. 

From elementary beam bending theory one finds that the shear 
force Q is given by 


Q = -Bw£’ 


(3) 


Since the shear strain is equal to w^, the shear force must also 
obey 


Q = GAw^ , 


#For the sake of simplicity, new equation numbering sequences 
are begun in each appendix. All equation numbers mentioned in 
this appendix refer to equations presented in this appendix. 
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density p 2 $ n (u)) of the fluctuating pressure due to turbulence on 
an EBF is given by Eqs. (5) and (10), but with the latter reduced 
by a factor of 4.4 (or 10 log 4.4 * 6.5 dB). 

Correlations . — If one approximates the pressure field as 
spatially homogeneous, as was done for the momentum annihilation 
case, one may write the wavenumber-frequency spectrum of the 
boundary-layer-like fluctuating pressure field as 

$ p (k,m) = p 2 $ n (to)$ x (k x )$> y (k y ) (17) 


where k, k x , k y have the same meanings as previously. The forms 
of the component spectra $ x and $ y and the magnitudes of the 
parameters that enter them may be determined on the basis of data 
given in Ref. 8 relating to the spatial correlations of the fluc- 
tuating pressures. One finds that 


Mk ) = 


L /tt 
x 


1 + (k - k ) 2 L 5 
x h 


L /tt 

* (k ) « 2 — 

y y 


1 + k L 


2t 2 


y y 


(18) 


(19) 


Here k^ denotes the hydrodynamic (or "convective") wavenumber, 
i.e., the wavenumber corresponding to pressure fluctuations 
that pass an observation point at the flow convection velocity 


k h = «/V c . (20) 


This convection velocity is related to the jet exit velocity U 
as 0 


V c * 0.45 U Q (21) 
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$ (k . 10 ) = 

P -* 

4> p (x» T ) = 


h/// 

111 


4> p (x» T ) 

$ (k .(jo) 
P 


-i(k • X-0.T) dKdydT 

(11) 

i(k -k-a.T) dk du 

x y 

(12) 


The pressure correlation <j>p(s,x) given by Eqs . (6) and (3) 

is isotropic; — i.e. , it is a function of only the separation 

s = [(xj - x 2 ) 2 + (y : - y 2 ) 2 ] 1/2 rather than of the coordinates 

(Xjjyj) and (x 2 ,y 2 ) of the two observation points. This isotropy 

implies that the wavenumber-frequency spectrum is a function of 

only the magnitude k = [k 2 + k 2 ] 1/2 of the wavenumber, rather 

than of the vector k; that is. 


$ ( k _ to ) = $ (k,w) 
P - P 


(13) 


Since the pressure correlation <(>p(s,T) consists of a pro- 
duct of a spatial and a time function, the wavenumber-frequency 
spectrum consists of the product of a frequency and a wavenumber 
function ; 


$ (k,w) = p 2 $ n ((ju)$(k) 


where here 


*00 = 


(2tt) : 


e -s/L e -ik -x dxdy> 


(14) 

(15) 


With the above relation for the separation s, 

o 

i f 


t(k) = 


2 TT 


" s/L J„(ks)ds - [1 + k 2 L 2 ] 3/2 


(16) 


where J is the zero-order Bessel function of the first kind, 
o 
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suggest that this correlation depends only on the separation 
distance s between the two observation locations, and not on 
the direction from one point to the other. Furthermore, these 
data indicate that the correlation length (i.e., the distance 
within which the correlation decays to a small fraction of its 
maximum value) is considerably smaller than the distances over 
which there occur significant variations in the mean velocity 
U or in the mean-square fluctuating velocity u T . Thus, one is 
justified in considering u to be a nearly spatially homogeneous 
field. Taking account of the exponentially decaying character 
of the correlation indicated by the data, one therefore may 
write 


u u 
1 2 


(t) = 


4> u (s,t) = 


u 


r -s/l -M/Tf 


(3) 


where L and T f denote the correlation length and the correlation 
time, respectively, of the fluctuating velocity. 

The data also suggest that one may estimate the two above- 
mentioned correlation parameters from 


L » 2.5 x 10 2 X (4) 

T f - 0.1X/U o (5) 


where U represents the jet exit velocity, 
o 

In terms of the previously indicated approximations, one 
may express the pressure correlation (Eq. 2) on the EBF as 

4) (x) = 6 (s, t) = 4p 2 U 2 <j> (s,x) (6) 

PjP^ P u 

For distances X of interest for realistic flaps (i.e., X/D * 9)» 
the local velocity U is very nearly equal to the exit velocity 
U Q or U c (see Figs. 14b and 15). 

Mean-square pressure. — From Eq . (6) one finds that the 

mean-square fluctuating pressure p 7 is given by 

^ = <J> (0,0) = 4p 2 U 2 u r » ( 4ql ) 2 (7) 
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APPENDIX A 


) 


MODELS AND ESTIMATES OF AEROACOUSTIC LOADS 
ON EXTERNALLY BLOWN FLAPS 


Jet Efflux Configuration 

A typical idealized jet efflux configuration is sketched in 
Fig. 2, showing a converging conical "potential core" surrounded 
by a spreading conical expansion region. It is expected that 
externally blown flap (EBF) structures will be inserted in this 
spreading region, in order to deflect the flow downward, so as 
to provide lift. The problem to be considered here consists of 
characterizing the fluctuating pressures that will act on the 
flap surfaces, so that one may estimate the corresponding vibra- 
tory stresses for design purposes. 

At first glance it may seem inappropriate to consider the 
flow produced by a fan-jet engine like that from a simple jet. 
However, one may conclude from Fig. 14 that the simple jet gives 
a reasonable approximation to the flow profile in those regions 
which EBF surfaces are likely to be. Whereas Fig. 14a indicates 
that a pronounced low-velocity fan-flow annulus may be discerned 
near the exit plane, Fig. 14b shows that the fan flow plays only 
a relatively minor role at locations several diameters from the 
nozzle plane. Indeed, if one plots the fan— jet data of Fig. 14b 
in non— dimensionalized form and compares it with similar data 
for an ideal circular jet (Fig. 15), one finds that the fan-jet 
profiles are quite consistent with profiles* corresponding to 
an ideal jet issuing from the core nozzle (with diameter D c ) 
with the core exhaust velocity U . 

If one assumes that the aeroacoustic noise field in the 
engine exhaust stream is, like the flow profile, similar to 
that for an ideal circular jet, then one may use the rather 
extensive data available for ideal jets to estimate the fluc- 
tuating pressures acting on an EBF inserted in the exhaust 
stream. 


*Such dimensionless mean velocity profiles are known to 
vary only little with the jet Mach number M. 
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CONCLUDING REMARKS; RECOMMENDATIONS 


The approach suggested here for estimation of fluctuating 
pressures associated with engine exhausts is based on extrapola- 
tion of nondimensionalized fluctuating pressure data obtained 
from simple jets and on interpretation of similar velocity dis- 
tribution data from measurements in the exhausts of a very lim- 
ited number of fan-jet engines. Clearly, the availability of 
data on the fluctuating pressure distribution in the exhaust 
of the engine to be used in any particular application may be 
expected to improve the characterization of these pressures and 
to" increase the confidence one has in fatigue life estimates 
based on these pressures. Comparison of predictions based on 
the approach suggested in this report with corresponding full- 
scale fan-jet engine fluctuating-pressure data would also serve 
as a useful check on the validity of the suggested approach. 

The response, stress, and fatigue life estimation approaches 
presented in this report follow the earlier literature in assum- 
ing only the fundamental mode of the structure to be of impor- 
tance. Although this assumption may lead to conservative designs 
and life estimates in many cases, one can easily visualize prac- 
tical situations where higher modes predominate. Such cases are 
particularly likely to occur with engine exhaust excitation, 
where the excitation pressures are correlated over small areas, 
have spectral peaks at frequencies considerably higher than the 
fundamental structural resonance, and convect along the struc- 
tural surface. Indeed, there also exists some experimental evi- 
dence that shows that higher structural modes play important 
roles in responses to flow excitation. Of course, the importance 
of higher modes in determining fatigue life is also enhanced by 
the higher fatigue damage accumulation rates associated with 
their higher resonance frequencies. Thus, it appears advisable 
to use the response, stress, and fatigue estimation approaches 
suggested here with some caution. Reexamination of these ap- 
proaches, and their extension to include appropriate higher mode 
responses, is recommended. 

Any but the most grossly empirical fatigue life prediction 
method must be based on information concerning how the .number of 
loading cycles that a structure can withstand varies with the 
fluctuating stress. The method suggested in this report is based 
on sonic fatigue data derived from tests on panel specimens of 
only one material for each panel type. In particular, the data 
pertaining to the fatigue of honeycomb-core sandwich structures 
is extremely limited. Thus, although one may expect the sug- 
gested prediction technique to yield good results for structures 
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Approximations for Small Angular Displacements , 
angular displacements, for which cos0 n = 1, sin©^ = 6 
above relations reduce to 

F xl = ”'1 X 1 + k x(l,0) l - X l ■ c " Be i ] + k x(l,2) [x l' x 2 +G + 
F yl - m ih + k y(l,0) [ h + B - c9 i ] + k y ( 1 ,2 ) ^- y 1 - y 2 

h i + F yi D " F xi E = h 0 i +k eu,o) [0 i “ 0 o ] + k e(i,2) [e i 

- k x(i s o) tx r c - B0 i ](B - ce i ) - k y(i,o) ty 

+ k x(l 2) [x l ■ x 2 + H +N +H0 1 + L0 2 ][H 
- k y(l,2)[ y l - y 2 - H - L + G6 1 +N6 2 ][G 

F x 2 = m 2 X 2 + k x(2,0) [x 2 ' K ■ J0 2 ] + k x ( 2 , 1) C x 2 " X 1 " G ' 
+ k x(2,3) [X 2- X 3 + R + U + Se 2 +T9 3 ] 

P y2 - m 2 y 2 + k y(2,O) [y 2 +J ' K0 2 ] + k y (2 , 1 ) [y 2 * y l + H ' 
+ k y(2,3) [y 2' y 3- 5 - T + R0 2 +l,0 3 ] 


- For small 

, the 

n ’ 

ii +110.^ + L©2 1 
L + G0 1 + N0 2 ] 

-0 2 ] 

1 +B-C0 1 ](C+B0 1 ) 

- G 0 x ] 

+ H0 1 ] 

- N - H0 . - L 2 ] 
f L - G 0 x - N 0 2 ] 
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m 2 X 2 + k x(2,0) C? 2 " ^0 ] (2,0) + k x ( 2 , 1 ) U 2 " ^1 ] ( 1 , 2 ) 


+ k , 0 

y (2 

F y 2 = m 2^2 + k x ( 2 ) 0) [n 2 _ n 0 ] ( 2 s 0 ) + k y (2 

+ k y (2 

M 2 + F y 2 P “ F x 2 Q = X 2®2 + k 0 ( 2 , 0 ) ( 6 2 " 9 0 

“ k x ( 2 , 0)^2 " V ( 2 , 0 ) 
+ k y ( 2 , 0) [n 2 " n 0 ] ( 2 , 0 ) 

- k x ( 2 , 1)^2 " ^ 1 - 1 ( 2 , 1 ) 
+ k y (2 , l )^ n 2 ” n i ' ] ( 2 , l ) 

_ k x ( 2 , 3) [C 2 " ? 3 ] ( 2 , 3 ) 
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,3) [? 2 " ^3 ] (2,3) 

,l) [n 2 " n l ] (l,2) 

,3) Cn 2 " n 3 ] (2 } 3) 

} + k 6(2,l) (0 2 " 6 1 ) 

+ k 0 ( 2 , 3) (6 2 " 9 3 ) 

[ Jcos0 2 - Ksin0 2 ] 
[-Kcos0 2 - Jsin0 2 ] 

[Lcos0 2 - Nsin0 2 ] 
[-Ncos0 2 - Lsin0 2 ] 

[-Scos0 2 + Rsin0 2 ] 


+k y(2,3) [T ’2' T1 3 ] (2,3) [Roos0 2 + Ssln0 2 ] 


whence 


^(n)(n,s) ^ (s ) (n ,s ) 


? (n) ” 5 (s) 


(n,s ) 


= x - x + X / xcos0 - X / \C0s6 

n s n(n,s) n s(n,s) s 


- Y , x sin9 + Y / x sinB 
n(n,s ) n s (n,s ) s 


( 100 ) 


(n)(n,s) n (s)(n,s) 


n (n ) ” n (s ) | (n,s ) 


= v — v + Y / x. cos 0 — Y / _ x cos 0 

^n ''s n(n,s) n s(n,s) s 


+ X , xsin0 - X , xSinO^ 
n(n,s) n s(n,s) s 


One may note that all k's, X’s, Y’s, as well as F x , F y , 

M, m, and I represent known or given quantities. On the other 
hand, the x's, y's and 0’s are the unknowns. Thus, one has 
three unknowns per flap element, as well as three equations 
of motion per element, so that one has as many equations of 
motion as there are elements. 


Equations for Illustrative Flap System 

Dimensions . — The rather formidable appearance of the 
foregoing equations is due to the somewhat intricate notation, 
which was introduced for the sake of generality. In order to 
obtain a clearer view of the nature of these equations (and 
of the meaning of the notation), it may be instructive to 
refer to the particular configuration sketched in Fig. 13. 

If the various capital letters indicated in that figure denote 
(positive) dimension values, then in terms of the previously 
introduced notation, 


H6 


displaced by the amounts ?,n in the x and y directions, respec- 
tively, if the flap center of gravity is displaced by the amounts 
x n and y n in the coordinate directions, and if the flap is 
rotated through an angle 0 n , where 


£ = x + XcosG - YsinO 
n n n 

(97) 

n = y + Ycos0 + Xsin0 
J n n n 


Spring Forces . — The difference between the displacements of 
the force-interconnection points on two flap elements, which 
points coincide when the flap system is in equilibrium, may thus 
be found from the difference between the £'s and p’s for the two 
points. This difference corresponds to the extension (or com- 
pression) of the interconnection springs, and therefore determines 
the spring forces. 


Dynamic Equilibrium of Flap Element . — Using the attachment- 
point coordinate designations* indicated in Fig. 11, one may find 
the following equations of motion for the nth flap element: 


F 

xn 

m x + 
n n 

k x (n ,n-l ) 

^(n) 


+ 

k x (n ,n+l ) 

? (n) 

F 

yn 

*• i 

my + 
rr n 

k y (n ,n-l) 

n (n) 


+ 

k y (n,n+l) 

n (n) 


? (n-l) 

? (n+l) 

n (n- 1) 
n (n+l) 


(n ,n-l) 


(n,n-l) 


(n,n+l) + 


(n,n+l) + **• 


(93a) 


*E.g., Y , , s denotes the Y coordinate, as measured in the 

& * n(n,n-l) 5 

system attached to the nth flap element, of the point at which 

the nth and (n-l)th elements are interconnected. 



Changes in the flap element's flexural rigidity and mass 
per unit length may be seen to have somewhat lesser effects 
on these two fatigue lives. On the other hand, the core 
density affects the core's fatigue life very significantly, 
with a 10 % increase extending the fatigue life by a factor 
of about 2.7. Again, the jet exit velocity is the most 
important jet parameter, with a 10 % decrease in U Q leading 
to increases in the skin and core fatigue lives by factors 
of 2.9 and 6.4, respectively. 
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If one again substitutes for $p from Eq. (8l), one obtains 


640 L : 


2.26 


7.05 \ / 

d 10.6 n 3.525 \ / 

/ U 3,525 p* 


05 y 10. 57 5 f 

oo 


L h 


10 10 t 7 ' 05 p d" 10 ! 6 

re f ref 


for f << f. 


3.1 


^ d 10 ' 6 (A/e) 7 - 05 f | 3,525 
l L 12 ' 1 ( y/B ) 1- 2 6 5 


\ I x 3,525 

1 l 7.05 y 17.6 2 5/ 

'Vo o 


for f >> f ip (94) 


Design Considerations. — Equations (90) and (93) show that 
the fatigue lives of facing sheets and honeycomb cores increase 
with decreasing exciting pressure spectral density (evaluated at 
the flap element’s fundamental natural frequency). This trend 
is as one would expect intuitively. Since the spectral density 
decreases with increasing frequency, as indicated by Eq. (9), one 
may obtain greater fatigue life by designing the flap element to 
have a higher fundamental resonance. In view of Eq . (55), a high 

fundamental resonance results from use of short unsupported spans 
L and of large stiffness/mass ratios B/y. 

As evident from Eqs . (92) and (94), reductions in L can re- 

sult in quite dramatic increases in fatigue life, provided that 
the flap element’s fundamental frequency f x is above the transi- 
tion frequency fip. If f x < f , then the facing sheet fatigue, 
life increase produced by a given amount of length reduction is 
somewhat less dramatic - and this length reduction may indeed be 
expected even to reduce the core’s fatigue life. 

From Eqs. (91) and (93) one may determine that 


^h 2 . 6*10 11 
If ~ B, 



n i r Q ref 

. c f 3, P (f - ) J 


1 , ■» 9 5 


(95) 
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a 2 . 1x10“ 4 , 


which indicates that stringers typically have much longer fatigue 
lives than. skin panels — in agreement with experimental observa- 
tions . 



Thus, one generally should first design the panels so that 
they have adequate fatigue lives, and then verify that the stringer 
design selected (usually on the basis of other than fatigue consid- 
erations) has a fatigue life that is no less than that of the pan- 
els . 


Honeycomb-Core Sandwich Flaps 

Facing Sheets (Skin) . — If one substitutes into the fatigue 
life expression of Eq. CT 7) the number of cycles to failure as 
given by Eq . (71) and the natural frequency expression of Eq . (55) 

— with a; = 2.36, to account for. boundaries that are neither sim- 
ply supported nor fully clamped, — and if one also uses the sec- 
ond stress expression of Eq. (68), which applies for the same 
boundary conditions, one finds that 



(90) 


In order to display the dimensional correctness of this ex- 
pression clearly, one may define an effective flap density 
Pf = p/A, an effective radius of gyr ation r = /B/EA, and an ef- 
fective longitudinal wavespeed c-p = /E / p -p . With these substitu- 
tions, one may rewrite Eq . (90) as 
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for f << f,p 
(87) 


for f >> f r , 


Design Considerations . — As is evident from Eqs . (80) and 

(86), the fatigue lives of skin panels and of stringers increase 
as the excitation — represented by the pressure spectral density 
$p(f) — decreases. This behavior is as one would expect intui- 
tively, of course. Since this spectral density decreases with 
increasing frequency, as indicated in Eq. (7), one should design 
the skin panel to have as high a fundamental resonance frequency 
f as possible. In view of Eq . (79), this implies that one should 
choose the largest admissible panel thickness h and the smallest 
panel edge length a. One might also consider choosing materials 
with large longitudinal wavespeeds cl, but most acceptable struc- 
tural materials have wavespeeds that differ by no more than about 
10% from each other, so that one stands to gain little by choos- 
ing alternate materials on this basis. 

Table I, which has been developed on the basis of the high- 
frequency part of Eq . (83), shows by what factors the fatigue life 

of a blown flap panel may be expected to change as the result of 
changing the various parameters. Thus, for example, one finds 
that a change in the skin thickness h by a factor of 1.50 (i.e., 
a 50% increase) would increase the fatigue life by a factor of 
71; similarly, decreasing the panel edge length a by 20% ’would 
lengthen the fatigue life by a factor of 1*1, whereas doubling the 
damping n would increase that life by a factor of 4.9. 

As evident from both Eq. (83) and Table I, relatively small 
changes in h and a can lead to quite considerable changes in pan- 
el fatigue life; the effects of changes in the other structural 
parameters are much less significant. Small changes in the dis- ■ 
tance X of the flap from the engine exit have relatively little 
effect on the fatigue life, and changes in the gas density (asso- 
ciated with exhaust temperature changes) that can occur with a 
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If the panel resonance frequency f is in a range where 
$p(f) varies little with changes in f, then Eq. (80) exhibits 
all of the dependences of the fatigue life on the panel para- 
meters. On the other hand, if <J> p (f) varies significantly with 
f, then the dependence of f on the panel parameters gives rise 
to additional effects. From Eqs, (8), (15) and ( 16 ) one finds 
that 

0.0115 Pq U qX for f << f T 

(81) 

0.29 p 2 U*/Xf 2 for f >> f T 
where the transition frequency f,p obeys* 

f T * 1.5 U Q /X. (82) 


Substitution of Eq. (8l) into ( 80 ) and use of Eq . (79) 

results in 



where the first parentheses enclose all relevant panel para- 
meters, whereas the second enclose the jet parameter terms. 

Stringers. — Again for the purpose of exhibiting the sa- 
lient parametric effects most simply, it is useful to consider 
the common case where the stringer length is the same as the 
greater of the two panel edge lengths, and where the spacing 
between stringers is equal to the shorter panel edge length. 
With b s = b and a s = a, assuming b/a > 3, and using Eq. (79), 
one may approximate Eq . (48) by 


*For the typical values of U Q = 750 ft/sec and X - 10 ft, 
one finds fip * 110 Hz. 
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Estimation Procedure . - In order to estimate the fatigue life 
of a honeycomb core sandwich flap element, the following procedure 
is recommended; 

1. Calculate the fundamental beam resonance frequency f x 
from Eq. (55). If the beam end conditions are not 
well defined, take a x = 2 . 36 . 

2. Determine the spectral density $n( f i) of the pressure 
acting on the flap at frequency f l from corresponding 
data, or estimate it from Eqs . (15) or (16). 

3. Find the maximum root-mean-square stress in the skin from 

Eq. (65), if the end conditions (and the associated fund- 
amental mode shapes) are well defined. Otherwise, find 
that stress by use of Eq . (68). Take p - 0.0^4, unless 

better data are available. 

4. Estimate the maximum root-mean-square shear stress in 
the core as the basis of Eq. ( 69 ). 

5. Calculate N f from Eq . (72) or (73) and find the skin 

fatigue life Tp from Eq . (7^0 for the confidence limit 

of interest. For materials other than 5052-H39 aluminum, 
multiply Tft by k m from Appendix C. 

6. Calculate from Eq. (75) or ( 76 ) and find the honey- 
comb core fatigue life T h from Eq. (77). For materials 
other than 7075-T6 aluminum, multiply Tf by k m from 
Appendix C. 

7. Take the effective fatigue life of the entire structure 
as the lesser of the values of Tp and T^. 
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Fatigue Life 


Facing Sheets. - From Fig. 51 of Ref. 9, which summarizes the 
results of a regression analysis of data obtained on panels with 
7075-T6 aluminum alloy facing sheets, one may deduce the following 
relation between the maximum skin stress a x and the number of 
cycles Nf that the skin can withstand without failing: 


log N f = -4.06 log I — 


or 



The discussion that follows Eq. (46) applies here again; the 
fatigue life T fC of the facing sheet corresponding to the -C% 

confidence limit may be found from 

TfC ' N fC /f . (71 ° 

where, of course, f represents the fundamental natural frequency 
of the sandwich beam and is found from Eq. (70) or (71) • 
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rms stress associated with beam vibrations in the fundamental 
mode as 
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1 


Ec ( u 


it 

rms 


max 


ECU 


rms 


<b" 

max 


( 64 ) 


where the primes indicate differentiation with respect to the 
lengthwise coordinate x, E denotes the modulus of elasticity of 
the skin material, and " represents the maximum absolute value 
of <r(x). max 

If one combines Eqs . (60), (6l), (62), and (64), one finds 
that the maximum root-mean-square skin stress obeys 


a 
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it J \ n ' 


(65) 


The similarity of the term under the square-root sign to the 
square-root terms appearing in Eqs. (29) and (47) is obvious. 

The term J <J>" depends only on the mode shape, and thus 
i max 

only on the boundary conditions. Since the mode shape, normalized 
in accordance with Eq. (57), is given (Refs. 23, 25) by 


4> ( x ) 


/2 sin(fTx/L) 


:osh(Bx/L) - cos(|3x/L) 


for simply supported ends 

( 66 ) 

y[sinh(Bx/L) - sin(Bx/L)] 
for clamped ends 


where 3 ~ 4.730, 


4tt/L 


Y ~ 0.9825, one may determine that 
12.6/L for simply supported ends 
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37.4/L 


for clamped ends . 


(67) 
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where L denotes the beam's length, B its bending stiffness*, and 
y its mass per unit length. The parameter a : represents a con- 
stant that depends on the boundary conditions; for a beam that is 
simply supported on both ends, a 1 = tt/ 2 * 1.57, and for a clamped- 
clamped (or free-free) beam, a 2 ~ 3-56. Since in a realistic flap 
element the boundary conditions are likely to be somewhere be- 
tween simply suppo rted and fully clamped, one may reasonably take 
* /( 1.57) ( 3'- 56 ) = 2.36 as a first estimate, in absence of better 
information . 

Resonant Response of Fundamental Mode . — In order to analyze 
the response of a uniform beam in its fundamental mode in general 
terms, it is convenient to introduce the mode shape <J)(x) associ- 
ated with that mode. One may then express any time-dependent beam 
deflection u(x,t) in the first mode (Ref. 33) as 

u(x,t) = U(t)<|>(x) • (56) 

The mode shape 4> ( x ) is defined physically only within a multipli- 
cative constant; of the various normalizations possible, the one 
chosen here (to facilitate use of available tables and references) 
is 


f L 

[ <f> 2 ( x ) dx = L . ( 57 ) 

'O 


It is well known that the dynamic response U(t) of any struc- 
tural mode is like that of a simple spring-mass-dashpot system 
with a mass equal to the modal mass 


M 


1 


r L 

y4> 2 (x)dx = yL 

Jo 


(58) 


exposed to a force that is equal to the modal force 


*For a homogeneous beam of a material with Young's modulus E and 
with a section having a moment of inertia I, the bending stiff- 
ness is B = El. 
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The following procedure is recommended for estimation of 
stringer fatigue life: 

1. Calculate the panel fundamental resonance frequency f 
and the pressure spectral density $p(f) as in the panel fatigue 
life estimation procedure. 

2. Evaluate the approximate maximum rms stress in the 

stringer from Eq . (48), using n = 10~ 2 , unless better damping 

data are available. 

3. Find the corrected rms stress estimate from Eq. (49). 

4. Calculate N r from Eq . (50) or (52) and find T n 

SO s u 

from Eq . (54). 

5. For materials other than 7075-T6 aluminum, multiply 

T _ by k from Appendix C. 
sC J m ^ 
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between the rivet lines of adjacent stringers, H the stringer 
depth (see Pig. 7)> and I an effective moment of inertia of 
the stringer cross-section, given by 


I 


'xx 


J xz 

X zz 


(48a) 


Here I xx denotes the centroidal moment of inertia of the stringer 
cross section about an axis parallel to the panel surface (see 
Fig. 8), I zz represents a similar moment of inertia about an axis 
normal to the panel surface, and I X z denotes a similar mixed mo- 
ment of inertia.* As previously, f denotes the fundamental reso- 
nance frequency of a skin panel, q represents its loss factor, 
and $p(f) denotes the spectral density of the fluctuating excita- 
tion pressure (at the frequency f). 

Correction of rms stress estimate on basis of test data . — 
Since the various assumptions involved in the derivation of Eq . 
(48) may represent rather poor approximations of conditions oc- 
curring in practical structures, one would expect predictions 
made on the basis of Eq . (48) to deviate from corresponding ex- 

perimental results. Comparison of such predictions with experi- 
mental data (Ref. 7) indicates that the experimentally observed 
root-mean-square stress a e on the average is related to the cor- 
responding calculated from Eq . (48) as 
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a 

e 


psi 



psi 



(49) 


The above relation was derived on the basis of calculated 
values ranging from about 150 to 3000 psi; its applicability to 
values outside this range remains to be established. 

Cycles to failure; survival probability; fatigue life . — 
Data presented in Pig. 44 of Ref. T indicates that the number N s 


•That is, I xx => / A z’dA , I zz = I A x 2 dA , I xz = / fl xzdA , 

where A represents the area of the stringer cross-section. Ap- 
pendix I of Ref. 7 gives expressions for these moments of inter- 
tia for zee, channel, and hat sections. 
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Although the time-variation of a randomly varing parameter 
like the panel stress is not a simple sinusoid, and one can not 
speak of cycles in the strictest terms, one may expect the panel 
vibrations to be approximately sinusoidal in time as long as they 
are dominated by a single mode — as was previously assumed in the 
response analysis. One may then consider the stress signal be- 
tween successive zero crossings as a half cycle, with the signal 
varying approximately like a sinusoid at the natural frequency of 
the system. Knowing this frequency and the number of cycles N 
that produce failure, one may calculate the fatigue life. 

For failure probability distributions that are symmetric 
about the mean, (50-C/2)$ of the samples fail under fewer stress 
cycles than the number corresponding to the -C$ confidence limit 
(Ref. 21 ). Thus, for example, if N = 10 7 cycles corresponds to a 
confidence limit of - 50 $ for a given panel design exposed to a 
given excitation, one may expect 25$ of all panels to fail at 
less than 10 7 cycles (i.e., one may expect 75$ of all panels to 
survive after 10 7 cycles). 

Fatigue life . — The fatigue life of a structure obviously 
must be defined in terms of a failure probability or similar 
statistical measure. Here it is convenient to use the fatigue 
life corresponding to the -C$ confidence limit, which one may 
find from 

T c = N c /f , (47) 

where N- is obtained from Eqs. (45) and (46) for the confidence 
limit in question and f denotes the natural frequency of the panel 
under consideration. 

In order to estimate the fatigue life of a given panel, the 
following procedure is recommended: 

1. Calculate the panel fundamental resonance frequency from 
Eq . (35), using the correction of Eq. (37) for curved panels. 

2. Determine the spectral density $p ( f ) of the pressure 
acting on the panel at the resonance frequency f from correspond- 
ing data or from Eq . (15) or (16). 

3. Find the maximum root-mean-square stress from Eq. (4l), 
modifying the result according to Eq . (42) for curved panels. 

Take n = 10" 2 , unless better data are available. 
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Fig. 5.2 of Ref. 8); the estimates tend to be too high - by as 
much as a factor of 5 at low stress values, and generally by less 
at higher stresses. Means for improving the estimates are not 
available at present. 

Below are tabulated values f R /f and c r R /a calculated from 

Eqs . (37) and (42), taking for b 2 /hR the greatest reasonable 

value. This value is b 2 /hR = A 2 (a/h) (a/R) ~ 175 A 2 , where 

rrid.A 1 1 icx a 

(a/h) = 500 is the greatest likely practical value for this 

max 

ratio, and where (a/R) = 0.35 is the greatest value of which 

ma x 

Eqs. (37) and (42) hold. For smaller values of a/h and a/R - 
that is, for smaller b 2 /hR - both f R /f and a R /o are nearer to 
unity. 

b/a 0.3 1.0 3.0 

f R /f 1.55 8.45 13.1 

a„/a 0.76 0.33 0.83 

n 

Panel loss factors . — Because the resoonses of panels to 
random excitation are dominated by the responses of resonant 
modes, the damping of a skin panel — as characterized by the loss 
factor n — is important in establishing the magnitude of its re- 
sponse and the associated oscillatory stresses. As evident from 
Eq . (29) and from the relations derived from that equation, the 
root-mean-square stress varies inversely as /rf. 

References 18 and 19 suggest a method for estimating the 
loss factors of panels with riveted edges, taking account of such 
parameters as rivet spacing, width of contact area, and panel 
wavelength (as a function of frequency). However, this method 
may be somewhat too cumbersome for preliminary design purposes. 
Abundant experimental evidence* indicates that for conventional 
aircraft structures (i.e., for structures not provided with 
special damping treatments), q differs little from 10 -2 . This 
value may therefore be taken as a reasonable estimate, unless 
measured data for the particular structure under consideration 
are available. 


*E.g., see Ref. 20. Reference 8 suggests n = 0.0085 for typical 
aircraft structures, based on values between 0.008 and 0.009 
reported for fuselage panels in Ref. 13, and on 0.0085 reported 
for tailplane panels in Ref. 14. The loss factors of the test 
panels investigated in Ref. 9 ranged between 0.005 and 0.009. 
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On the other hand, if one uses the expression (35) 
to more realistic boundary conditions, one finds 


a = 


1.10 


V 


an 




corresponding 


(40) 


In Ref. 9 there is derived yet another expression for a, on the 
basis of a regression analysis of experimental data for 
1 £ b/a <_ 3 and 190 £ b/h <_ 500 , arranged in nondimensional 
grouns of variables deduced from ideal clamoed-edge panel analy- 
sis. This expression (when rewritten in consistent units) cor- 
responds to 


a = 


0.24 


V 


ClV 77 

an 



, - 0.0 6 


p- 1.6 8 


(41) 


It is instructive to compare the stress estimates one ob- 
tains from the three foregoing equations. Clearly, Eqs. (39) 
and (40) differ only in the functions of the aspect ratio b/a 
they involve and (slightly) in their numerical coefficients. If 
one evaluates the aforementioned functions (Fig. 8), one finds 
that for practical values of the aspect ratio the two functions 
differ by about ten percent in magnitude and exhibit very nearly 
the same trend. Because of the differences in the magnitudes of 
these functions and in the coefficients of Eqs. (39) and (40), 
the stress one calculates by use of Eq . (39) is higher by a fac- 
tor of about 1.2 than the stress estimate one finds from Eq, (40). 


The function F 1 ' 68 , which appears on the right-hand side of 
Ea. (4l) obviously depends more strongly on b/a than the func- 
tion F^ 2 associated with Eq . (39); this stronger dependence is 

evidenced by the slightly steeper slope of the F 1 ' 68 curve in Fig. 

3 / 

8. The function F 1 ' 68 also is found to exceed F /2 by about 21% 

for b/a = 1, and by about 67%> for b/a = 10. However, Eq. (41) 

differs from Eqs. (39) and (40) also in the exponent on (b/h) 

and in the added n -0 ' 06 term. For the typical value of n = 10 -2 , 

the latter term amounts to 1.32. For values of (b/h) between 

100 and 500, corresponding to the panels used in tests which 
served as the basis of development of Eq. (4l), the value of 

(b/h)/ 4 lies between 3.17 and 4.72. Thus, on the average, 

(b/h) 7/4 n- 0 * 06 F " 1 - 68 /(b/h) 3/2 F" 3/2 - (4.0) (1.3) (1.4) — 1 ~ 3.7. If 
one multiplies the numerical coefficient 0.24 by this value, 
one obtains 0.90; comparison of this result with the coefficient 
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and that of a clamped panel obeys (Ref. 9) 

f ~ il JR F(b/a) 
1 c 3 T m ab 


( 30 ) 


(3D 


where m denotes the mass per unit area of the panel and R its 
flexural rigidity, and 

F(b/a) = [2 + 3b 2 /a 2 + 3a 2 /b 2 ] 1/2 . (32) 

For homogeneous panels 


D 


Eh 3 

12(1 - v 2 ) 


( 33) 


m = p g h , ( 3^ ) 

where E represents Young's modulus, v Poisson's ratio, and p 3 
the density of the structural material. One may thus estimate 
the fundamental natural frequency of a skin-stringer panel from 


f = Jr f = — - G(b/a) 
V s c 2 


( 35) 


where c^ = /rVp s represents the longitudinal wave velocity* in 
the panel material and 



/ 2 \ 

Va ( 

G (b/a) = g- 

! I 1 ♦ fr) p(b/a) 

« ) 
'1 


.25 


for b/a = 1 
for b/a >> 1 


(36) 


*For most structural metals, one may take c^ = 2 * 10 5 in/sec 
with adequate accuracy. Since v 2 << 1 typically, Poisson's Ratio 
does not appear in these approximate expressions. 
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be well justified In many practical cases. Once the panel motions 
have been determined, the sonic fatigue analysis method proceeds 
to determine the stresses Induced In the stringers as the result 
of the loads Imposed on them by the vibrating panels, and it then 
compares these calculated stresses with experimental stress and 
fatigue data. 


Skin 


Relation between dynamic and static stress . — In Ref. 6 
analysis of the response of elastic structures to random pres- 
sure fields is discussed in general terms, and simplified re- 
sults are presented for the case where: 

(1) One mode predominates in the frequency range of 
interest . 

(2) The excitation pressure is in phase over the 
entire structure of interest. 

(3) The spectrum of the excitation does not change 
rapidly in the vicinity of the resonance fre- 
quency of the dominant mode. 

It is shown in Ref. 6 that If the foregoing conditions hold, 
then the root-mean-square stress o induced at a given location 
in a structure (or panel) by a random pressure field may be ex- 
pressed in terms of the stress a 0 induced at that same loca- 
tion by a uniformly distributed static pressure of unit magni- 
tude as 


a 


V 


TT 

2n 


f $ (f ) 
n p n 


o 


(29) 


Here n represents the structural loss factor (of the dominant 
mode at its resonance), f n denotes the (cyclic) resonance fre- 
quency of the dominant mode, and $p(fn) represents the spectral 
density of the exciting pressure at the frequency f n . 

Resonance frequencies of flat rectangular panels . - A meth- 
od for determining the natural frequencies of multi-bay systems, 
taking into account the flexural and torsional stiffnesses of 
the various stringers, is presented in Ref. 7. An alternate 
method, applicable to structures with many equi-spaced identical 
stringers between flexurally stiff frames, is summarized in Ref. 
8. These methods, however, are relatively complex; — perhaps 
too complex for preliminary design purposes. It is likely also 
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corresponds to the previously chosen reference values and to the 
reference frequency f „ = 100 Hz. Substitution of Eq . (20) 

then yields re 



T + 460 

-20 log ° 52Q (28) 
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Overall fluctuating pressure level . — The overall fluctua- 
ting pressure level Lqa is a logarithmic measure ( ex pressed in 
decibels) of the mean-square fluctuating pressure p 2 , defined 
as * 


L qa = 10 log (p 2 /p^ ef ), dB (17) 

where p re p is a reference pressure. Substitution of Eq . (13) 
into the above then permits one to write 

L 0A = c + 20 log (p o /P ref ) + 40 log (V U ref ) <l8> 


where P re f and ^ref are reference values of the density and 
velocity , and 

C = 20 log (0-24p ref Uj >ef /p ref ) > dB • (19) 

If one chooses P re f = 0.0735 lb/ft 3 (corresponding to air at 
at room temperature and one atmosphere), U r ef = 750 ft/sec (a 
typical core engine exhaust velocity), and Pref = 2.9 x 10 _9 psi 
(= 0.0002 microbar, the international standard reference value 
for acoustic pressures), one obtains C * 177.5 dB. 

Since the density of a gas is inversely proportional to 
its absolute temperature, one may replace the density ratio of 
Eq. (18) by a corresponding temperature ratio and write 

/ U \ T + 460 

l oa “ 177 • 5 + 2,0 log ur^l " 20 log — 9 2 0 — » (20) 


where T q represents the temperature of the exhaust stream in °F. 


Pressure spectrum level . — The spectrum level of the 
fluctuating pressure is a logarithmic measure (expressed in 
decibels) of the spectral density $ (f ) , defined as 


L g = 10 log 


_^(f) 


^Pjref^" 


( 21 ) 


where <J> p >re f(f) is an appropriate reference value, usually taken 
as p 2 ref /Hz. 


*A11 logarithms in this report are base 10. 
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V p 2 /< 


which show how \ p 2 /q varies in the stream-wise and cross-wise 
directions along a plate inc line d at various angles to an im- 

finds that Vp 2 /q does not exceed 0.1 for 


flow-normal 

smaller 


ringing j et , one 

X/D _> 9 and for angles between the plate surface and 
greater than 30°. This value of 0.1 is considerably 

than the value of Vp 2 / q = 111 * 11(0.12) = 0.H8 one obtains from 
Eq . (5) for a normally impinging jet with the near-maximum tur- 
bulence intensity I = 0.12. 


The frequency-spectral density $p(ai) of the fluctuating 
pressures associated with nearly tangential flows may again be 
approximated by Eqs. (7) and (8). The pressure cross-correla- 
tion function here is more complex, however, being character- 
ized (see Appendix A) by different correlation lengths (or 
"eddy decay scales") L x and Ly in the stream-wise and trans- 
verse directions, with 


L * 13.5 V /w 
L a 2.0 V /w 

y c 

where 

V » 0 . H5 U 
c o 


( 11 ) 


( 12 ) 


represents the convection velocity of the flow. 


Design Pressures and Pressure Levels 

Maximum mean-square pressure . — Although one may use the 
data shown' in Pigs. 3-6, together with the previously given 
equations, to estimate the fluctuating pressures that occur at 
any specific location, one usually need not consider all this 
detail for design purposes. By inspection of Figs. 3 and ^ 
one finds that for 8 £ X/D < 20, corresponding to typical lo- 
cations where EBF surfaces may be expected to be placed normal 
to the flow, the turbulence intensity does not exceed 0.12. 
Since one also may note that in the high-turbulence region 
(i.e. for r/D <_ 1 ) the turbulence intensity decreases slowly 
with X/D, approximately according to Eq. (3), one may choose 
I = 0.12 for general conservative design purposes. 
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varies along the axis of a jet, and Fig. 4 shows how this in- 
tensity varies along the radial coordinate. Here U 0 represents 
the jet exit velocity (which, for a f an-j et engine is taken to 
be the core engine exit velocity) and u 2 denotes the mean-square 
axial fluctuating velocity. 

As is evident from Fig. 3, the intensity I on the jet axis 
is at approximately its maximum value of 0.11 at X/D ~ 10. Fig- 
ure 4 shows that for X/D _> 8 , I does not exceed approximately 
0.12. From examination of the peak values of Fig. 4 one may de- 
termine, in fact, that for X/D > 9, the maximum value of I obeys 

I * 0.165 - 0.0044 X/D (3) 

max 

and occurs at a radial coordinate r p ea ^> given by 


P g ak * 0.15 | - 1.0 (4) 

The velocity and pressure fluctuations within the potential 
core typically are much smaller than those in the flow outside 
the core. Thus, for the regions of interest with respect to 
blown flans, I * 0.12 may be expected to represent an upper 
bound suitable for conservative design purposes. 


Pressure Fluctuations on Flap Surfaces 

Normally Impinging Jets . — For flap surfaces on which the 
jet flow impinges essentially normally, one may take the momen- 
tum flux in the flow to be annihilated at the structural su rf ace 

With this assumption, the mean-square fluctuating pressure p 2 is 

found to be related to the mean-square fluctuating velocity u 2 
(see Appendix A) as 


p 2 = 4p 2 U 2 u 2 * ( 4ql ) 2 , (5) 

where p denotes the local fluid density and U the local mean ve- 
locity. For most locations of interest for EBF ' s , the local ve- 
locity U is nearly equal to the exit velocity U and the local 
gas density p differs little from the density p° at the exit. 

With these assumptions one obtains the above indicated approximate 
equality, where 
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VARIATION OF SURFACE PRESSURE FLUCTUATION ON A FLAT PLATE DUE TO 
NORMAL JET IMPINGEMENT (0 = 0°) WITH X/D, AT SEVERAL LOCATIONS ALONG 



AXIAL TURBULENCE INTENSITY, I=</j2/u ( 



FIG. 4. VARIATION OF AXIAL TURBULENCE INTENSITY WITH DISTANCE 
FROM JET CENTERLINE (FROM REF. 4). 
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Response Calculation 


Beam . — For the purpose of illustrating the general technique, 
it is again convenient to discuss the response of a uniform beam. 
Treatment of more complex structures involves no different concepts 
only more intricate calculations. 

From the definitions of the admittance and of the pressure 
spectral density, one finds that the mean-square velocity of the 
beam at any given frequency oj x obeys 


v 2 ( 0 ) 


i 


) 





) | Y ( k , a) x ) | 2 


dk 


( 8 ) 


Thus, if one knows both 4>(k,u) and Y(k,oi), one may calculate the 
mean-square velocity at any frequency. 

Because the two functions appearing in the integrand typically 
have pronounced peaks, these peak values usually dominate the re- 
sponse, permitting one readily to obtain simple approximations for 
the integral. Figure 22 illustrates this behavior for a one-dimen 
sional structure subject to simultaneous jet noise and turbulent 
flow excitation. 

The upper part of this figure shows a sketch of the pressure 
spectral density function, which for this simultaneous excitation 
consists of the sum of the functions due to the separate excita- 
tions (see Figs. 19 and 20). Thus, the density function here ex- 
hibits one peak at the acoustic wavenumber k D (at the frequency 
uj ! under consideration), and a second peak at the hydrodynamic 
wavenumber k^ . Also shown is the square of the admittance func- 
tion (see Fig. 21) and the product of these two functions, which 
is the integrand of Eq. (8). The relative magnitudes and loca- 
tions of these peaks, of course, depend on the relative strengths 
of the two excitation sources, on the structural damping, on the 
structural mass and stiffness, on the flow speed relative to the 
speed of sound, and also on the frequency. 

A k,w plot like that shown in the lower portion of Fig. 22 
permits one to obtain considerable insight into the variation of 
the positions of these peaks. As previously discussed in connec- 
tion with Figs. 19 and 20, the pressure spectral peaks occur at 
the acoustic and hydrodynamic wavenumbers, which correspond to 
intersections of the acoustic and flow velocity lines, respective- 
ly, with the fixed frequency line w = w . Similarly, the admit- 
tance peak occurs at the wavenumber k^ = k b (w 1 ); this wavenumber 
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Surface Pressures due to Jet Noise . - The typical behavior 
of the spectral density associated with surface pressures pro- 
duced by jet noise is indicated in Fig. 20. One may note that 
the contour plot here is much like that of Fig. 19, except that 
the contour lines here cluster about the to = kc G line, indicat- 
ing that most of the energy associated with jet-noise-induced 
surface pressures travels at the speed of sound. 

If all of the sound from the jet would impinge on the sur- 
face at grazing incidence, all of the contour lines would col- 
lapse upon the c 0 line. In practice, however, a flap surface 
is subject to sound arriving with a distribution of angles of 
incidence, as well as to acoustic nearfield components, so that 
there occurs a distribution of energy about the sound speed 
line. Nevertheless, the peak in the spectral density at a giv- 
en frequency oij occurs at the acoustic wavenumber k 0 = ui /c 0 . 

The peak of the $(k,w) hill is a noteworthly feature. It 
represents a concentration of fluctuating-pressure energy that 
generally is important for structural fatigue and noise consid- 
erations, and it occurs at the jet wavenumber k-- * 2Tr/Dj , where 
Dj represents the jet diameter. 

Spectral Characterization of Structural Response 

Admittance of Beam . — It is instructive to illustrate ap- 
plication of the spectral response characterization for a one- 
dimensional system, such as a beam or one-dimensional plate de- 
forming in flexure; generalization to two-dimensional systems 
then can be accomplished relatively simply. The well-known 
equation of flexural motion of a uniform beam is 

B + m ^2- = p ( x , t ) , (4) 

3X 1 * 3t 2 

where y represents the beam's lateral displacement, B denotes 
its flexural rigidity, and m its mass per unit area. On the 
right-hand side there appears the exciting load per unit length 
p(x,t),which is a function of the axial coordinate x and time t. 
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the most extreme case, the structural response calculation then 
involves the same mathematical process as for more conventional 
space-time descriptions — namely, convolution of a spatial struc- 
tural kernel with a load cross-correlation function. At worst, 
the STSD approach leads to no more complexity than other approaches. 


Spectral Description of Random Pressures 

Spectral Densities . — From classical Fourier analysis it 
follows that a pressure wave that is harmonic over all space 
and time can be represented (in terms of the usual complex vari- 
able notation) by an amplitude and an exponential phase factor. 
It is also well known that a general pressure function can be 
represented as an infinite sum of such harmonic waves; this rep- 
resentation is called the Fourier transform. 


A pressure wave that is periodic in space and time, that 
travels with velocity U in the positive x-direction , and that 
passes any fixed point at the (radian) frequency w, thus may be 
described by 


where p(k,w) denotes the amplitude of the wave and 

k = w/U ( 1 ) 


is known as the wavenumber. Generalized harmonic analysis per- 
mits any arbitrary random pressure p(x,t) to be represented by 


p( x,t ) 



v i(kx-wt) 
w)e 


dkdw , 


( 2 ) 


where p(k,w), the density of the amplitude distribution in k and 
a) space, is the Fourier transform of p(x,t). 

From the square of the magnitude of this density one may 
determine the so-called pressure spectral density 
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APPENDIX D 


AN INTRODUCTION TO THE SPACE-TIME SPECTRAL DENSITY APPROACH 
TO EXCITATION AND RESPONSE CHARACTERIZATION 


Introducti on 

As has been pointed out in the main body of this report, 
the responses of flap structures to fluctuating pressures de- 
pend not only on the temporal (frequency) characteristics of 
these pressures, but also on their spatial characteristics. 

The presently available sonic fatigue analysis and design ap- 
proaches, including those summarized in the main body of this 
report, avoid the complexities associated with accounting for 
the spatial characteristics by making the assumption that the 
exciting pressures are uniformly distributed over the structure 
under consideration and that the structure’s most significant 
response occurs in its fundamental mode. Although this assump- 
tion often leads to conservative designs, one may readily show 
that it need not in all cases — and, indeed, there exists some 
flight data (e.g.. Refs. 15 and 16 ) that indicate that the 
structural responses are not described adequately by the sonic 
fatigue analysis approaches in current use. 

This appendix serves as a brief introduction to an approach 
which should be able to provide a logical framework for the guid- 
ance of data acquisition programs for load and response charac- 
terization, and which also may be expected to lead to more real- 
istic predictions. This approach, which makes use of space-time 
spectral density (STSD) concepts, has been developed quite ex- 
tensively for dealing with the vibrations of ship structures 
induced by sound and flow (and with the underwater sound radi- 
ated by these vibrations), and has been applied to such problems 
with considerable success. As discussed below, it has the addi- 
tional advantages of being relatively simple in concept, of lead- 
ing to little computational difficulty, of permitting great ex- 
perimental simplification and savings, and of leading to increased 
accuracy . 

Conceptual and Computational Simplicity . - For spatially 
homogeneous excitation, the STSD approach permits one to calcu- 
late the structural response as an integral over the product of 
two functions, of which one characterizes the fluctuating pres- 
sure, the other, the dynamic stiffness of the structure. The 
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given by Eq . (6). The exponent (3, of course, depends on the 
component — i.e., on which of the aforementioned equations 
apply. The table at the end of this appendix gives the correc- 
tion factors corresponding to these equations and to a wide 
variety of aluminum alloys. 

Estimation of Fatigue Stress Ratio . — Unfortunately, S-N 
curves are available for only a few alloys. For alloys for 
which no fatigue data are available, one may use the rough 
approximation that 



( 7 ) 


where Y represents the yield stress of the material (Ref. 3b). 
As evident from the table at the end of this appendix, this 
approximation is very close for some materials, but may be 
about 20$ too high or too low for others. Nevertheless, in 
absence of better information, one can do no better than to 
use the above relation. The correction factors given in the 
table below are based on fatigue stress ratios where these 
are available, and on yield stress ratios otherwise. 

Materials Other Than Aluminum . — It should be noted that 
the procedure suggested here for aluminum alloys cannot readily 
be extended to other materials, unless their S-N curves have 
the same slopes (on a log-log plot) as those for aluminum. 
Unfortunately, most other materials have different slopes and 
many — notably steels — have segments of greatly differing 
slopes. For such materials, further analysis and/or experi- 
mental investigation is required. 


67 



APPENDIX C 


FATIGUE LIFE CORRECTIONS FOR ALUMINUM ALLOYS 


The various fatigue life estimates presented in the main 
text were obtained on the basis of experimental data on struc- 
tures made of only one kind of material — namely, 7075-T6 
aluminum for skin/stringer structures and for the facing sheets 
of honeycomb sandwich structures, and 5052-H39 aluminum for 
honeycomb cores. This appendix suggests how one may correct 
the estimates pertaining to the aforementioned materials so 
as to obtain corresponding estimates for other aluminum alloys. 

Similarity of S-N Curves . — The fatigue behavior of mate- 
rials generally is described by so-called "S-N" curves, which 
are plots of the fully reversed stress amplitude S versus the 
number of stress cycles N at which a specimen fails when sub- 
jected to cyclic stress of that amplitude. When plotted on 
log-log scales, the S-N curves for most aluminum alloys appear 
very nearly like parallel straight lines, at least in the low 
stress and large N region [E.g., see Ref. 33 and Table 3*3*l(c) 
of Ref. 34] . Although the classical S-N curves are obtained 
from experiments where the stress amplitude is held constant 
(for each data point), whereas the S-N curves represented by 
Eqs . (45), (52), (71) and (75) correspond to random stress 
variations with a given mean-square value, one may expect the 
latter log-log curves for various alloys to be parallel, if 
the former are parallel. 

If one assumes that the root-mean-square stress a that 
different alloys can withstand for a given number of cycles 
is proportional to the "fatigue stress" S of the material, then 
for two different materials (indicated by subscripts 1 and 2), 



For materials that exhibit a definite endurance limit (i.e., a 
stress amplitude that the material can withstand essentially 
for an unlimited number of cycles), one would be inclined to 


*For the sake of simplicity, new equation numbering sequences 
are begun in each appendix. All equation numbers mentioned in 
this appendix refer to equations presented in this appendix. 


65 



where G denotes the effective shear modulus of the core (consid- 
ered as a continuum) and A represents the beam's cross-sectional 
area. By equating Eqs. (3) and (4) and integrating the result 
one may determine that 


w 


s 


-B „ 
AG W b 


(5) 


Fundamental Natural Frequency . — Substitution of Eq. (5) 
into Eq’. (2) yields - * 

llT !lir _ y , ,2 ( B It \ 

W b B \ b “ AG b/ • 

For a deflection given by 


TTX / 7 \ 

w b = w 0 sm — , U) 


which corresponds to the fundamental mode of a simply supported 
beam, Eq. (6) yields 


where 


03 





tt 2 B \ ^ 
L 2 AG / 



( 8 ) 


(9) 


may be recognized as the classical expression for the natural fre- 
quency, corresponding to the case where shear effects are neglected 
(i.e., where AG, the shear stiffness per unit length, is assumed 
infinite ) . 

Flexural and Shear Stresses . — The flexural stress in the 
outer-most skin fiber, taken to be a distance c from the neutral 
axis, may be found from elementary beam theory to be given by 

a = -Ecw b = Ec(tt/L) 2 w fl sin(iTx/L) , 


( 10 ) 



The parameters L x and Ly represent the eddy decay scales in the 
flow and transverse directions, respectively, and obey 


13.5 V c /M 

(22) 

2 V c /M 

(23) 


The pressure field described by Eqs. (17) to (23) is the 
same as that associated with a turbulent boundary layer, except 
for the numerical constants that appear in Eqs. (21) to (23). 
Detailed interpretations of such fields are given in Ref. 29; 
some salient features are summarized below. 


Because the hydrodynamic wavenumber appears in the de- 
nominator of Eq . (l8), the wavenumber spectrum $ x (k x ) is asym- 

metric in k x . This asymmetry represents a mean convection of 
the fluctuating pressure field in the positive x-direction with 
a velocity V c . The fact that L x and Ly have finite values ac- 
counts for the decay of the correlations witn increasing separa- 
tion. Inverse variation of L x and L with frequency (in keeping 
with the "similarity hypothesis" of nef. 32 ) accounts for the 
decrease in the correlation lengths with increasing frequency. 

Figure 18 indicates qualitatively how $p(k,w) varies with k x 
and ky at a fixed value of w. The peak value of $p(k,w) occurs 
at k x = k h , ky = 0. Because L x > Ly , the shape of the $ p (k,w) 
function is much more elongated in the ky than in the k x direc- 
tion. Since the peak occurs at the nonzero wavenumber k^ , there 
is a region centered around the wavenumber vector k = (kh,0) over 
which $p(k,w) has relatively high values. This region, which is 
important in relation to vibration response, is known as the "con- 
vective region" of the wavenumber plane; its location and extent 
clearly are frequency-dependent. 
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Comparison of estimated and measured spectral density , - The 

frequency spectral density p 2 $ n (w) of the fluctuating pressure 
measured (by means of a small flush-mounted microphone) on a half- 
scale EBF model , at a point near the axis of a cold jet (at X/D =9, 
with Up ~ 750 fps, pressure ratio 1.4) is shown in Fig.l6, together 
with an estimate based on Eqs . (5), (7) s ( 10 ) and on Fig. 3. Con- 
sidering the courseness of the various approximations that underlie 
this estimate, the agreement between it and the data is quite 
reasonab le . 


Estimation on Basis of Boundary Layer Flow 


Turbulent boundary layer pressures . — The foregoing estima- 
tion approach was based on the assumption of momentum annihila- 
tion, and thus in essence assumed flow impingement essentially 
normal to the flap surface. However, since at least at some lo- 
cations on the flap the flow is essentially parallel to the sur- 
face (see Fig. 17), it is useful to consider a model of the fluc- 
tuating pressure on an EBF that resembles that for a turbulent 
boundary layer. 


Extensive data on pressure fluctuations produced at the 
surface of a flat plate by a jet impinging at various angles 
are reported in Ref. 5. Figures 5 and 6 reproduce some of this 
non-di me n s ionali zed data, showing how the root-mean-square pres- 
sure Vp^varies along the plate surface. One may observe that 
for fixed X/D and 0 the mean-square pressure p 2 varies slowly 
with the distance x along the plate, indicating that the fluctu- 
ating pressure field is spatially inhomogeneous to a slight ex- 
tent. As evident from Fig. 6, this inhomogeneity decreases at 
increasing distances X from the jet exit plane. 


may find that 
layer pressures 


From Figs . 5 and 6 one 
an upper bound for boundary 
figurations (9 ~ 30 ° > X/D = 9) 
hand, one finds that for normally impinging 

^Vq ~ 0.44. 


VpVq » 0.1 provides 
typical EBF con- 


on 


From Eq. (8), on the other 

flow (with I = 0,11) 


Examination of data of Ref. 5 pertaining to the frequency 
spectra of fluctuating pressures on the surface of the test plate 
(see Fig. 5 for test geometry) indicates that the normalized fre- 
quency spectrum $> n (u>) of the fluctuating pressures for the turbu- 
lent boundary layer case has the same shape as that for the nor- 
mal impingement case. Thus, one may estimate that the spectral 
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where 


* = I pU o 


(8a) 


denotes the dynamic pressure at the jet exit and 



(8b) 


represents the axial component of the turbulence intensity on 
the axis at the axial distance X of interest. Figure 3 shows 
how I varies with X; note that for the typical location X/D = 9 
for an EBF, one finds I » 0.11. 


Normalized correlations and frequency spectrum . — From 
Eqs . ( 1) and ( 6 ) one finds that the normalize d time-correlations 

of pressure and velocity obey 


♦ (0/0 <t> u ^,T) 

4> p (0,0) = <J> u (0,0) = e 


The normalized frequency spectrum $ n (u) is the Fourier transform 
of the normalized correlation, so that the normalized frequency 
spectrum of pressure (or velocity) obeys* 


$ ( co) 
n 



T f /1T 
1 + m z T z 


(10) 


Wavenumber- frequency spectrum . — The general space-time 
correlation <f>p(x,t) is closely related to the corresponding 
wavenumber-frequency spectrum <J>p(k,a)), where x = (x,y) is a 
two-dimensional position vector (representing the Cartesian co- 
ordinates x and y of the observation point in the plane of the 
flap and k = (k x ,k y ) is the corresponding wavenumber vector. 

For SDatially homogeneous and temporally stationary fields 
(Ref/ 29), 


*Unless otherwise 
taken from negative to 


indicated, all integrations 
positive infinity. 


are to be 
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Estimation on Basis of Momentum Flux 


Pressure-momentum relation . — At any given location in the 
jet, the axial momentum flux' may be written as p(U + u) 2 in terms 
of the fluid density p, the mean axial velocity U and the fluctu- 
ating axial velocity u at the location of interest. If this mo- 
mentum flux is entirely annihilated as the jet impinges on the 
EBF surface, then the fluctuating component p of the pressure 
acting on the surface is of the order of the fluctuating compo- 
nent of the momentum flux. Thus, for the usual case where the 
fluctuating velocity component is much smaller than the steady 
component (u << U), one obtains 


p a 2 p uU 


(l)* 


Cross-correlations. — The cross-correlation <}> v ^ (t) of the 

P i P 2 . 

pressures p x and p 2 at two different locations on the flap, de- 
fined as 

_ (t) = < p (t)p (t + t) > 

p i p 2 

thus is related to the corresponding velocity cross-correlation 
<{> (t) = < u i (t)u 2 (t + t) > 


as 

♦„ „ CO * u CO . (2) 

P j P 2 12 

In the above expressions, the brackets <...> denote averaging 
with respect to time t; t represents a time interval, and U x 
and U 2 denote the mean velocities at the two locations of 
interest . 

In order to obtain some simple estimates readily, it is 
convenient to consider locations near the jet axis (i.e., near 
point A of Fig. 2 ) in a plane normal to that axis. Data for the 
velocity correlation <b (0) for this special case (Refs. 27,28) 

U , U o 


*For the sake of simplicity, new equation numbering sequences 
are begun in each appendix. All equation numbers mentioned in 
this appendix refer to equations presented in this appendix. 
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that are closely related to those whose fatigue data were used 
in development of the technique, its reliability is much reduced 
for other structures. Accumulation of a more extensive data 
base is recommended, particularly for the types of configura- 
tions and materials likely to be used for future externally 
blown flaps. 

Because of the great potential utility of the space-time 
spectral density approach, it is recommended that application of 
this approach to EBF and related problems be pursued vigorously, 
both in relation to characterization of the fluctuating pres- 
sures produced by the impingement of engine exhaust on flap sur- 
faces and in relation to the estimation of structural response 
spectra. 
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'1(0,1) 

= -C 

X 2(0,2) 

= -K 

X 3(2,3) 

1(0,1) 

= B 

Y 2(0,2) 

= J 

Y 3(2,3) 

1(1,2) 

= G 

X 2(l,2) 

= -N 

X 0 
a3 

1(1,2) 

= D 

Y 2 ( 1 , 2 ) 

= L 

Y 0 
a3 

al 

= D 

X 2 ( 2 , 3 ) 

= R 


al 

= E 

Y 2 ( 2 , 3 ) 

= -S 




X a2 

= P 




Y a2 

= Q 



Equations of Motion . — The equations of motion for this 
configuration then may be written as 

F xl = m l x l +k x(l,0)^l ■ 5 ohl,0) + k x(l,2)^l - e 2-'(l,2) 

F yl = B l»l +k y(l,0) tn r n 0 1 (l,0) + k y(l,2) [r1 l ~ n 2 ] (l > 2) 

M 1 + F yl D - F xl E * h S l + k e(l,o) [0 l - S 0 ] + k e(l,2) [0 l- e 2 ] 

- k x(l,0) [5 l * V(l,O) [Boos0 l - Csln6 l ] 

- k x(l,2)[q - { 2 ] (l,2) [ - HcoslS l +Gsin0 l ] 

+ k y(l ) O) [n l- ,1 O ] a,O) [ - Ccos0 l- Bslne i ] 

+ k y(l,2) [r, l -’ 1 2 ] (l,2) CGc0se i + Hslne i ] 
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M+FX -FY =ie+k fi , , v ,v - u . 

n yn an xn an n n 0(n,n-l) V n n-1 


( 6 n 6 n-l) 

+ k 0(n,n+l) ( 9 n ' 6 n+l) + * * 


- k 


x (n ,n-l ) 


^(n) “ ^(n-1) 


, s Y , , N cos0 +X , , N sin0 

(n,n-l) n(n,n-l) n n(n,n-l) n 


- k 


x (n ,n+ 1 ) 


C (n) " ? (n+l) 


/ . n s Y / ,-,nCos0 +X / vSin0 

(n,n+l) n(n,n+l) n n(n,n+l) n 


+ k 


y(n ,n-l ) 


n (n) ~ n (n-l) 


, n n X , , v cos0 -Y , , N sin0 

(n,n-l) n(n,n-l) n n(n,n-l) u n 


+ k 


y (n ,n+l ) 


n (n) n (n+l) 


, , n , X , ,, A cos0 -Y , , sSine 

(n,n+l; n(n,n+l) n n(n,n+l) n 


(98b) 


Here m n represents the mass of the nth flap element, and I n 
denotes the moment of Inertia of that element about its center 
of gravity . F xn and F yn denote the unsteady components of the 
aerodynamic force acting on the nth element in the x and y co- 
ordinate directions, respectively, and M represents the aero- 
dynamic moment. Furthermore k x ( n s ) represents the stiffness 
of the spring, in the x direction^ connecting the nth and sth 
elements; k y(n s) represents the stiffness of the corresponding 
spring acting in the y direction; kQ( n s ) denotes the (rota- 
tional) stiffness of the rotational spring connecting these 
two elements. Also, if one lets 5( n )(n s ) and n (n)(n s) re- 
present the x and y displacement components of that attachment 
point between the nth and sth elements which is located in the 
nth element, one finds from Eqs.(97) that one may write 


^(n) (n,s ) 
n (n) (n,s ) 


x + X / N cos0 - Y , N sin0 
n n(n,s ) n n(n,s ) n 


y + Y , ^cos0 + X , >.sin0 
J n n(n,s) n n(n,s) n 


( 99 ) 
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GENERAL TWO-DIMENSIONAL EQUATIONS OF MOTION 
OF MULTI-ELEMENT FLAP SYSTEM 


Ideal i z a t i o n 

Masses and Springs . — In order to idealize a multi-element 
flap system as a two-dimensional (planar) dynamic system, it is 
convenient to consider each flap element (consisting of airfoil, 
guide-rails, and supporting structures) as a rigid body, and to 
take each such rigid body as connected to each other body at a 
single point by one set of springs — each set consisting of a 
single spring acting in the x-direction, one acting in the y- 
direction, and one acting rotationally . Thus, the flap system, 
as modeled in two dimensions, reduces to an array of planar 
rigid bodies interconnected by springs, as shown in Fig. 10. 

A schematic representation of the interconnecting springs is 
shown in the lower part of that figure . 

One may obtain the spring constants (ane Lytic&lly or experi- 
mentally) and the effective attachment point locations by dis- 
connecting all attachments except the one of interest, holding 
all bodies fixed except the one of interest, applying a force or 
moment, and observing the resulting displacements. 

Coordinates Attached to Flap Elements; Notation . —It is 
convenient to select a Cartesian coordinate system attached to 
each flap element, with the origin of this system located at 
the element’s center of gravity, and with the system’s X-axis 
aligned parallel to the x-axis of an inertial reference system 
(attached to the wing) when the flap is in its static equilibrium 
configuration. This coordinate system serves to locate the 
various force-application points (i.e., the aerodynamic force 
locations and the interaction spring attachment points) on the 
element with respect to the element's center of gravity. Figure 
11 shows the coordinate system on a typical flap element and 
indicates the notation used in the present analysis. 


Equations of Motion 

Displacement of General Point on Flap Element . — One may 
readily find (see Fig. 12 ) that a typical point P, whose coor- 
dinates are (X,Y) in the system attached to the nth flap, is 
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and from Eqs. (55) and ( 8 1 ) one may determine that for f >> f , 


n i r a ref 

c fV f i> 


4.8 



( 96 ) 


For typical orders of magnitude for the various parameters, one 
finds that 



so that, taking the middle value of B from Eq . (73), 


h ~ 2 . -6 xi o^ 1 ^ p5) 2.99 (50 ) 4 ' 06 (1) 10 -6 \ 7,05 [ ( 4 . 8)10 -2 (46 )] 49 5 

2 f 3.4xio 9 \10 3 / 

= 4. 6xio -9 . 

Thus, one generally would expect core shear fatigue failures to 
occur long before facing sheet failures.* In designing a flap 
element it thus appears logical first to select a core that 
has the required fatigue life, and then to verify that the 
facing sheet will endure at least for the same time span. 

Tables III and IV, which have been derived from the parts 
of Eqs. (92) and (94) that pertain to flap elements with high 
natural frequencies, indicate the factors by which the facing- 
sheet and honeycomb-core fatigue lives change as the result of 
changes in the various structural and jet parameters. Clearly, 
the one most significant structural parameter is the unsupported 
span length L; a mere 10$ decrease in L may be expected to 
increase the fatigue life of the skin by a factor of about 
4.5, and that of the core by about 3.6. 


*No comparable experimental data appear to be available. 
Such data as are available (Ref. 9) pertain to panels, rather 
than beams, and are affected by stress raisers (e.g., fasteners) 
that reduce the skin fatigue life. 
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Equations (90) and (91) exhibit all structural parametric 
dependences for the low-frequency case, where $p(f x ) is essen- 
tially independent of frequency. But, since ^(fj) does vary 
significantly at the higher frequencies, the aforementioned equa- 
tions need to be modified. If one substitutes for $ p (f) from 
Eq . (8l) and again uses Eq. (55) with = 2.36, one obtains 



As previously, the first set of parentheses in each expression 
encloses the structural parameters, the second the jet parameters. 


Core. — By substitution of the cycles-to-failure relation 
of Eq"! (76) and the natural frequency express Jon of ¥ x . (55), 
again using otj = 2.36, into the fatigue life equation, Eq . (77), 

and by using also the stress expressions of Eqs . (68) and (69), 

one may find that the honeycomb core fatigue life obeys 



The second form of this expression involves the parameters intro- 
duced after Eq. (90) and is presented here in order to demonstrate 
the dimensional consistency of this result. 
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given engine also are likely to have only minor significance. 

On the other hand, the exit velocity U Q is of great importance; 
a mere 10 % increase in this velocity can reduce the panel fatigue 
life to about one-third of its original value. 

Table II is analogous to Table I, but pertains to stringers. 
From Table II, which is based on the part of Eq. ( 87 ) that ap- 
plies for f >> fij, one may observe, for example, that doubling 
of the panel edge length a (which is also equal to the spacing 
between the stringers) increases the stringer fatigue life by 
a factor of 1.91, whereas doubling the stringer length b reduces 
the fatigue life to 0.47 of its former value. 

By comparing Eqs . ( 83 ) and ( 87 ) or Tables 1 and II one finus 
that the panel fatigue life is much more sensitive to parameter 
changes than is the stringer fatigue life. It is also evid« nt 
that increases in h and cj,, as well as decreases in a, serve to 
increase the panel fatigue life, whije they result in reductions 
in the stringer fatigue life. 

From Eqs. ( 80 ) and (86) one ,ay determine that the ratio of 
oanel to stringer fatigue life oocys 



If one introduces $p(f) as given by Eq. (8l) for f >> fry , if one 
substitutes for f from Eq . (79), and if one takes I = H 3 h/2 
(which corresponds to an I-beam of height and flange ’width H, 
with flange thickness h), one finds that one may approximate the 
above expression by 



for the purpose of making order-of-magnitude estimates. Substi- 
tution of the middle values of B and B s given in Eqs. (44) and 
(53) and of typical orders of magnitude of the various ratios 
then leads to 
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S 0.014 


o 
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Hb 2 
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V 


heJ^TT 


n 


(84) 


The "theoretical" maximum stress given by Eq . (48), how- 
ever, must be corrected according to Eq . (49), if one desires 

a better representation of the actual (experimentally observed) 
maximum stress o e . Substitution of Eq. (84) into (49) indicates 
that that stress obeys 


a 


e 


a 


re f 


1.5 


Hb 2 

|he L yn 

1 0 

re f 

V '' j 


(8b) 


where a _ = 10 3 ds 1, as before, 
re f - J 

By combining Eqs . (85), (79), (52) and (54) one obtains 
the stringer fatigue life T s as 



= 0. 48 B. 


<ho L ) 


a 2 / 1 ° re f ' 

1-274,1.096 \ H 


0.548 



2 7 4 


(86 ) 


where the first form again groups the parameters to display the 
dimensionless correctness and the second shows the parametric 
dependences more clearly. 

If one again uses Eq . (81) to account for the dependence 
of $ D (f) on the panel resonance frequency f, which frequency 
again may be approximated by Eq. (79), one may find that 
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DEPENDENCE OF FATIGUE LIFE ON JET AND STRUCTURAL PARAMETERS 


Skin-Stringer Flaps 


Skin Panels. — In order to display typical parametric de- 
pendence's convenient ly , it is useful to focus on commonly used 
skin-stringer configurations that have aspect ratios b/a > 2 . ‘3 
For such configurations, the function F(b/a) of Eq . (32) may 

be approximated by /3(b/a) and the function G(b/a) of Eq . (36) 
may be taken as equal to 0.69. 


If one takes the maximum panel stress to be given by Eq. 
(39), with the coefficient 1. 18 replaced by 0.90 in accordance 
with the discussion presented after Eq. (41), then one finds 
by use of the above indicated approximations that 


o 


0 . 395 


a 2 1 


f ) 


h -n 


V 2 


(78) 


and that the fundamental natural frequency cf the skin panel 
obeys 


f * 0.69 hc T /a 2 . (79) 

Substitution of Eq . (78) into (45), and substitution of the 

result and of Eq. (79) into Eq. (47) yields the following ex- 
pression for the panel fatigue life: 



The first form of this equation groups the parameters in a man 
ner that displays its dimensional correctness, whereas the sec 
ond form indicates the parametric dependences more directly. 



Core . - Only very few core shear fatigue failure data points 
appear to be available; the data do not suffice for regression 
analysis, and thus one can not establish confidence limits. 

Perhaps the best one can do at present is to accept the de- 
sign data indicated in Pig. 86 of Ref. 9, although the basis for 
that figure is not indicated. From the curves in that figure one 
may deduce the following relation between the number of cycles N, 
that will induce failure in 5052-H39 aluminum alloy honeycomb r: 
core, the maximum rms core shear stress t, and the core density d: 

log N h = 12.76 - 7.05 log^) + 10 . 5 8 log|g— j (75) 


or 


h 


5.7 


10 


12 ter te)‘ 


0 "6 


(76) 


Here t = 1 psi is a reference value of shear stress and d „ = 
rel ref 

1 lb/ft 3 is a reference value of density. 

Since not enough data are available for the determination of 
confidence limits, it appears reasonable to assume that the fore- 
going expressions pertain to the 0 % confidence limit. In order to 
estimate the numbers of cycles to failure corresponding to the 
-50% and - 95 % confidence limits (probably conservatively), one may 
mutliply the value of obtained from Eqs . (75) or (76) by 0.4 
and 0.1, respectively.* 

The honeycomb fatigue life T^q corresponding to the confi- 
dence limit C, of course, may be calculated from 


T 


hC 


N hC /f . 


(77) 


*These factors correspond approximately to the ratios of the con- 
stants in Eq . (53). The factors corresponding to Eqs. (46) and 
(73) are larger than those given here, hence would lead to less 
conservative estimates. 
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If one uses the geometric average (21.7/L) of the above values, 
in view of the uncertainty of the actual boundary conditions, one 
may rewrite Eq. (65) as 


a 


i 


0.69 


ecE ./ V 7 ^ 

UL 2 f^ V n i 


0.19 


ecE L / % (f 1 } 


( 68 ) 


where the latter expression was obtained by substitution of Eq . 

( 55 ) with a j = 2.36. 

Maximum Root-Mean-Square Shear Stress in Core . — In Appendix 
B it is shown by means of an analysis that parallels that pre- 
sented in Ref. 26 that, for a simply supported sandwich beam, the 
ratio of the maximum shear stress t in the core to the maximum 
tensile stress in the skin obeys 


ID t 

T _ ttB s 

a j ' LAcE 4 11 "L 


(69) 


where A represents the cross-sectional area of the beam. The 
approximate equality applies for a beam with a rectangular cross- 
section, with skin of thickness t g ; this approximate expression 
may suffice for the evaluation of a rough estimate in cases where 
not enough information is available to apply the more complete 
expression . 

For beams with other than simply supported boundaries, the 
simple analytical approach of Appendix B does not work and results 
like the above cannot be obtained readily. It is therefore sug- 
gested that Eq. (69) be used for estimation purposes, regardless 
of what the boundary conditions are. 

Loss Factors. — The available data pertaining to loss factors 
of honeycomb core sandwich structures are extremely limited. 
Reference 9 reports test results for about 30 different panels 
vibrating in their fundamental modes. Their loss factors were 
found to lie between about 0.03 and 0.05. and to be comparable to 
a value of about 0.04 measured on panels obtained from aircraft 
development programs. 

In absence of more directly applicable data, an estimate of 
r) = 0 . 04 appears to be reasonable. 
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p(x,t)<()(x)dx = ep Q (t) 


(59) 


Fj(t) 


J A 


rL 

<p(x)dx , 

Jo 


and with a modal stiffness 


kj = ( 2 tt f j ) 2 Mj 


( 60 ) 


In the above expressions, p(x,t) represents an arbitrary pressure 
distribution over the beam, p 0 (t) denotes a pressure that is 
spatially uniform, and e represents the beam width (see Fig. 9)- 
Note that here M is equal to the total mass of the beam. 

For a simple spring-mass-dashpot system subject to random 
excitation, one finds that the root-mean-souare displacement U 
obeys (see Refs. 23, 24) 


U 


rms 



( 61 ) 


where represents the spectral density of the force F 2 ( t ) and n x 
denotes the loss factor of the system, i.e., of the beam in its 
first mode. Since the spectral density of the force is pro- 
portional to the mean-square force, is related to the spectral 
density $ of the pressure p (t), in view of Eq . (59), as 

J P * - 0 J * 


<J> = e 2 J 2 $ 

F i p 


( 62 ) 


where 


J 
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rL 

(j) ( x) dx 
o 


(63) 


Maximum Root-Mean-Square Stress in Skin . — For a given amount 
of beam flexure, the greatest skin tensile and compressive stresses 
occur in those fibers that are farthest from the beam's neutral 
surface. If c denotes the distance from the neutral surface to 
the farthest fiber (see Fig. 9), then one may write the maximum 
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FATIGUE LIFE OF HONEYCOMB-CORE SANDWICH FLAPS 


Overview of Estimation Approach 

The fatigue of flat, rectangular panels of honeycomb sand- 
wich construction is discussed in Ref. 9 on the basis of classi- 
cal thin-plate theory. Earlier data cited in that report indi- 
cate that this theory yields good approximations to observed 
vibration and stress responses associated with the fundamental 
panel mode, and that shearing of the core plays no important role, 
unless this core is very flimsy. 

A flap element, however, may be expected to behave more like 
an end-supported beam than like an edge-supported panel. The 
analytical results available for panels thus do not apply to flap 
elements directly, although one may hope that honeycomb panel 
fatigue data will also be useful for honeycomb beam fatigue life 
estimation. The following paragraphs, therefore, first summarize a 
corresponding beam analysis and then apply related available fa- 
tigue data to develop a fatigue life estimation approach. 


Beam Response 

Resonance Frequency . — In order to simplify the analysis , 
it is useful as a first approximation to assume the fluctuating 
excitation pressure to be uniformly distributed over (one surface 
of) a flap element, and to consider only the response of the first 
mode of that element modeled as a uniform beam — in a similar man- 
ner somewhat analogous to that used in skin-stringer panel analyses 
or honeycomb sandwich panel analyses (Ref. 9). 

The resonance frequency of the first mode of a beam is given 
(e .g. , Ref. 22) by* 


f 


i 




(55) 


*This expression is based on the assumption that shear effects are 
negligible. As shown in Appendix B, the finite shear stiffness 
of a beam reduces its natural frequency, but in most practical 
cases this reduction is insignificant. 
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of cycles that a stringer survives before failing Is related to 
the experimentally observed maximum rms stress a e in it accord- 
ing to 


log N g = -2.74 logl 


ref 


+ log B 


(50) 


where o „ = 10 3 psi is a reference stress value as before, and 
ref K 


log B 


[6.98] 


f -95% ) 
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[7.94 J 

( 1 
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Equations (50) and (51) may also be rewritten as 


N = B ( a ,/a . „ ) 
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(52) 
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(53) 


The discussion and the relation between confidence limits 
and failure probability presented (in relation to panel failures) 
in the paragraphs following Eq . (46) apply equally well to string- 
er failures, as does the discussion of fatigue life. In analogy 
to Eq. (47), the stringer fatigue life T sC corresponding to the 
~C% confidence limit obeys 


T „ = N _/f 
sC sC 


(54) 


where N s q is found from Eq . (50) or (52) for the confidence limit 

of interest and f, it should be recalled, denotes the natural 
frequency of the panel. 
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from 


4. Calculate M„ from Eq. (45) or (43), and find T P 
Eq. (47). 

5. For materials other than 7075-T6 aluminum, multiply T P 
by k m from Appendix C. 

Stri ngers 

In typical skin~stringer structures, fatigue failures of 
stringers usually occur at the clip attachment (where the stringer 
is joined to the frame or bulkhead), because of the presence of 
stress raisers in that location. Because of the general complex- 
ity of the problem, little analytical work has been done on 
stringer fatigue, and since stringer failures generally occur in 
the interior of practical structures, there appears to exist no 
quantitative field data. Reference 7 contains the most defini- 
tive available analytical and experimental information; it is on 
that report that the following discussion is based. 

Analytical estimate of maximum root-mean •■square stress . — 

The analysis of stringer stresses presented in Ref. 7 is based on 
the following assumptions: (1) The total force acting on a stringer 

corresponds to the net shear force (integrated distribution minus 
corner reactions) that acts at the edge of a simply supported 
panel, which is deflecting in its first mode, in response to a 
pressure that is uniformly distributed over the panel, but vary- 
ing randomly in time. (2) The force acting on a stringer is dis- 
tributed uniformly along its length and acts on the rivet line . 

(3) The maximum stress in the stringer occurs in flexure at the 
clip attachment point, where the stringer is taken to be clamped 
with respect to bending. 

With these assumptions one finds that the maximum root-mean- 
square stress in a stringer obeys* 
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(48) 


where b denotes the stringer length (which is usually, but not 
necessarily always, the longer panel edge length;, a the distance 


*This relation follows from Eq . (67) of Ref. 7. However, there 

the numerical coefficient, which here is 2 3/ y3'n' 7 ^ ~ 0.0171, was 
erroneously omitted. 
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It should be noted, however, that air flow along a panel may 
extract energy from the panel vibrations, and thus Increase the 
effective structural loss factor — or that this flow may feed 
energy Into panel vibrations (under conditions approaching panel 
flutter), and thus decrease the effective loss factor. At pre- 
sent there is available no means for estimating this effect, and 
one can do little better than to evaluate it on the basis of 
experimental measurements. 

Cycles to failure; survival probability . — One may expect 
that the number of stress reversals a panel can withstand de- 
creases as the stress amplitude increases. Related test data, 
corresponding to skin-stringer panels of 7075-T6 aluminum alloy 
exposed to random noise, are given In Fig. 34 of Ref. 9, together 
with curves representing various statistical confidence limits. 
Later test data (Ref. 7) were found to fall within these same 
confidence limits; the design nomographs given In Refs. 7 and 9 
are based on these confidence limit curves. 

From the curves of the above-mentioned Fig. 34 of Ref. 9 one 
may determine that the number N of cycles that a panel survives 
before failing is related to the maximum root— mean-square stress 
a according to 


log N 5: -4 . 60 log 


(°ref ) 


+ log B , 


(43) 


where a „ = 10 3 psi is a reference stress value, and 
ref 
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Alternately, one may express the above relation as 


N = B(a/a ref )--‘ < * 60 


( 45 ) 


with 
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1. 18 of Eq . (39) indicates that the stress estimates one obtains 

by use of Eq . (4l) on the average are about 30$ lower than those 

one obtains on the basis of Eq. (39), and, similarly, to be about 
20$ lower than those one finds from Eq. (40). 


It is important to note, however, that Eq . (4l) was derived 
on the basis of experimental data (Ref. 9) on test panels with 
aspect ratios a/b between 1.0 and 3.0 only, so that the validity 
of this relation for larger aspect ratios remains uncertain. 
Furthermore, the test data points (see Fig. 69, p. 13B of Ref. 9) 
exhibit a good deal of scatter, with a large number of the points 
deviating considerably from the regression line. It thus is not 
clear whether the use of the somewhat more complex Eq. (4l) is 
justified instead of Eq. (39) with a reduced coefficient that 
makes this equation correspond more closely to the available data 
for b/a £ 3.0. 


In view of the fact that Eq . (4l) has gained some accep- 

tance, has been reduced to nomograph form, and has been com- 
pared with some experimental data (though net well docu- 
mented) other than that on the basis of which it was derived, 
it seems logical to retain it for stress estimation purposes. 
However, for the purpose of studying trends and parametric de- 
pendences, the simpler Eq. (39), with the coefficient 1.18 re- 
placed by 0.90, is likely to be advantageous. 


Maximum root-mean-square stress in curved panels . — On the 
basis of analytically developed expressions, in which empirically 
derived corrections have been included, the maximum root-mean- 
square stress op (at the middle of the straight edge) in a cy- 
lindrically curved panel with radius of curvature R has been 
found (Ref. 9,7) to be related to the corresponding stress in a 
similar flat panel as 
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(42) 


where, as before, A = b/a = length of curved edge/length of 
straight edge and f/f R is given by Eq. (37). 

The applicability of this relation is limited to 
0.3 £ b/a £ 3.0, a/h >, 100, and a/R £ 0.35. On the whole, 
stress estimates for curved panels obtained on the basis of Eq. 
(42) or corresponding nomographs correlate more poorly with 
test data than do similar estimates for flat panels (e.g., see 
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Resonance frequencies of cy lindrically curved panels . — The 
fundamental natural frequency f R of a cy lindrically curved panel 
with radius of curvature R is related to the natural frequency f 
of a flat panel with the same thickness h and edge dimensions 
as * 



0 . 006(b 2 /hR) 2 
A 4 + 0 . 6lA 2 + 1 






(37) 


where A = b/a, and a denotes the length of the flat edge of the 
panel, and b the length of its curved edge. The above relation 
was developed semi-empiri cally , on the basis of experimental 
data on structures with realistic boundary conditions, and is 
valid only for h/a £ 1/100 and for aspect ratios in the range 
0.3 < b/a < 3.0, and for a/R <. 0.35. 


Maximum root-mean-square stress in flat panels . — In a 
simply supported panel, the maximum flexural stress associated 
with uniform loading or with the first vibratory mode occurs at 
the panel center. In a rectangular panel that is clamped on all 
edges, the corresponding maximum stress occurs at the middle of 
the longer edge. In practical skin-stringer structures, panel 
fatigue failures typically occur along the edges, at the rivet 
line or at the ends of stringer flanges or doublers (Refs. 9, 
17); thus, the panel stress associated with fatigue corresponds 
more closely to the maximum stress in a clamped panel than to 
that in a simply supported panel. 


The maximum flexural stress c Q induced in a clamped panel 
with b £ a by a uniformly distributed static pressure of unit 
magnitude is given (Ref. 9) by 


°o max 



(38) 


where F = F(b/a) is given by Eq. (32). 


If one substitutes the foregoing for a Q into Eq . (29), and 
if one takes the natural frequency f to be equal to that for a 
clamped plate, one obtains the maximum rms stress as 


a 


1.185 


J5 


(39) 


*This expression results from Ref. 9, if a misprint in that 
report is corrected (see Ref. 12). Note: The expression appear- 

ing in Ref. 7 also is opviously in error. 
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that they give results whose precision is much greater than 
necessary, in view of the considerable uncertainty in (1) the 
estimated pressure spectra, (2) the validity in any practical 
case of the assumptions involved in the development of Eq. (29), 
and (3) the boundary conditions operative in practical struc- 
tures. Furthermore, fatigue data have been accumulated only 
for the lowest modes, so that application of this data to fa- 
tigue prediction for any given structure of a material or con- 
figuration different from those for which data is available, or 
to higher modes, is likely to introduce greater errors than 
those due to the use of simpler, less precise, resonance fre- 
quency estimates. 

Many measurements of the random responses of panels (e.g.. 
Refs. 6, 9—12) have shown these responses to be dominated by 
the fundamental panel mode. Data on realistic aircraft struc- 
tures (e.g., Refs. 13, 1*0 and related analyses (summarized in 
Ref. 8) have indicated that the responses of skin— stringer con- 
figurations generally* are dominated by modes in which each pan- 
el vibrates in a fundamental mode corresoonding to a boundary 
condition (at each edge) that lies between the fully clamped 
and the simnly supported. It is therefore reasonable to focus 
on the fundamental panel mode, and to omit the more complex 
higher modes from consideration. 

Reference 9 presents (on p. 224) a curve that summarizes 
the experimentally observed variation with aspect ratio of the 
fundamental resonance frequencies of rectangular panels of skin- 
stringer configurations. Although one may use this curve for 
estimation purposes, an analytic approximation to it will prove 
useful for determining how the fatigue life of the panel depends 
on the various parameters. Inspection of this curve, together 
with the corresponding curves for panels that are simply sup- 
ported and for panels that are clamped on all four edges, re- 
veals that the experimentally observed frequencies are very 
nearly equal to the geometric average of the resonance frequen- 
cies for the two ideal boundary condition cases. 

The fundamental resonance of a simply supported panel of 
thickness h and edge lengths a and b is given by 


*However, some data are available (Refs. 15, 16) which 
show that the fundamental panel mode response does not always 
predominate . 
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